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PREFACE 

This  report  describes  the  development  of  e  new  approach  to  modeling 
turbulent  flow.  Heat  transfer  in  a  fluid  is  also  addressed.  This  effort  is  in 
response  to  an  initiative  from  Klrtland  Air  Force  Base  calling  for  novel  methods 
of  modeling  incompressible  turbulent  flow  in  convectively  cooled  systems.  The 
approach  taken  has  been  to  derive  a  Hamiltonian  which,  preserves  all  of  the 
physics  contained  in  the  Navier-Stokes  equation,  and  to  diagonalize  this 
Hamiltonian  (i.e.,  transform  the  variables  to  variables  which  are  preserved  with 
respect  to  the  Hamiltonian).  This  Hamiltonian  can  be  used  as  the  starting  point 
for  a  number  of  different  calculations.  For  example,  it  can  be  used  to  obtain  a 
partition  function  that  will  be  used  to  describe  the  statistical  properties  of  a 
turbulent  system.  Indeed,  this  is  an  approach  favored  by  the  author.  The 
Hamiltonian  can  also  be  used  as  the  basis  of  the  newer  class  of  calculations 
referred  to  as  "Renormalization  Group"  and  "Chaos"  descriptions.  Hence,  it  may 
be  used  in  models  which  are  sometimes  referred  to  as  "deterministic"  as  well  as 
those  which  are  known  as  "probabilistic."  This  is  natural,  since  the 
Hamiltonian  itself  does  contain  the  same  physics  as  the  Navier-Stokes  equation. 

Rationale 

The  rationale  for  the  adopted  approach  is  that  it  has  worked  well  in  other 
areas  of  physics  which  share  with  turbulence  the  properties  of  nonlinearity  and 
a  large  number  of  degrees  of  freedom.  These  are  the  key  features  of  what  are 
known  as  many-body  systems.  Examples  of  systems  in  which  the  field  theory  or 
many-body  approach  has  proved  extremely  useful  are  nuclear  physics,  high-energy 
physics,  quantum-electrodynamics,  and  condensed  matter  physics.  In  all  of  these 


fields,  the  Hamiltonian  is  separated  into  the  sue  of  a  linear,  basis  part,  and  a 
nonlinear,  "perturbing”  part.  Solutions  for  the  full,  exact  system  are 
expressed  in  terms  of  integrals  involving  solutions  to  the  linearized  system. 
Usually  this  involves  some  form  of  perturbation  theory,  which  immediately  raises 
the  question  of  how  rapidly  a  perturbation  series  will  converge,  or  even  whether 
it  will  converge  at  all.  This  is  especially  true  of  turbulent  flow,  in  which 
closure  of  equations  involving  moment  expansions  of  the  velocity  components  is  a 
common  problem.  In  the  above  areas,  the  representation  of  the  complicated 
integrals  which  arise  in  perturbation  series  by  simple  diagrams  has  been  a  great: 
aid  in  summing  perturbation  series,  sometimes  to  infinite,  order.  This  has 
worked  even  in  cases  in  which  the  series  converges  very  slowly,  or  even 
diverges.  The  Hamiltonian  described  in  this  report  has  perturbations  which  go 
as  the  Mach  number,  and  so  perturbation  expansions  are  expected  to  converge  well 
for  subsonic  flow. 

The  usefulness  of  diagrams  to  represent  complicated  integrals  goes  beyond 
Its  great  simplification  of  a  bookkeeping  function.  The  diagrams  have  a 
psychological  appeal  that  aids  the  physical  intuition  of  the  theoretician  in  his 
formulation  of  a  particular  calculation  involving  interactions  between  various 
degrees  of  freedom  of  the  system  (e.g.,  fluid  elements,  or  normal  modes  of  the 
linearized  system).  In  the  above  fields,  the  various  degrees  of  freedom  of  the 
system  are  throught  of  as  particles,  and  interactions  are  described  in  terms  of 
scattering  events  between  particles.  For  example,  an  oscillation  of  a  sound 
field  at  a  particular  frequency  may  be  regarded  as  a  phonon.  This  viewpoint 
gives  a  deeper  meaning  to  the  term  many -body  system.  We  see  then  that  a 
many- body  description  of  a  problem  is  helpful  in  performing  diagram  calculations 
of  perturbation  expansions  and  highlights  parallels  with  other  systems  which  have 


been  extensively  analyzed.  Furthermore,  it  provides  a  means  of  making  a  generic 
analysis,  sinoe  calculations  are  made  in  terms  of  the  normal  modes  of  the 
linearized  system,  without  actually  specifying  explicitly  what  those  modes  are. 
The  boundary  conditions  and  explicit  solutions  are  included  at  the  end  of  a 
generic  calculation.  This  feature  may  facilitate  comparisons  of  the  behavior  of 
systems  with  different  geometries  and  flow  conditions. 

Overview 

The  various  aspects  of  treatment  of  a  system  based  on  field  theory  are 
shown  in  block  form  in  Figure  1.  Host  calculations  begin  with  the  full 
Hamiltonian,  H,  as  a  starting  point.  The  Hamiltonian  is  expressed  in  terms  of 
fields  which  are  the  solutions  to  the  partial  differential  equation  of  interest, 
in  this  case  the  Navier-Stokes  and  heat  budget  equations;  hence  the  name  field 
theory.  Note  that  at  the  outset  field  theory  gives  a  generic  description.  The 
fields  are  the  wavef unctions,  or  solutions  to  the  partial  differential  equation; 
calculations  are  done  in  terms  of  these  wevefunctions  without  explicitly  solving 
for  them.  Many  calculations  can  be  done  in  terms  of  a  Poisson  bracket  (P.B.) 
formulation.  Poisson  bracket  calculations  have  the  advantage  of  being 
independent  of  a  particular  coordinate  system,  and  are  very  powerful.  The  above 
mentioned  many-body  systems  are  all  quantum  systems,  and  the  analogous 
equations  to  P.B.  relations  in  those  systems  are  commutator  relations. 

It  is  emphasized  that  the  formulation  described  in  this  report  is  entirely 
classical.  At  no  point  is  a  quantum  description  used.  Although  the  formulation 
closely  parallels  field  theory  descriptions  of  quantum  systems,  such  as  the  use 
of  Poisson  brackets  instead  of  commutators,  the  integrity  of  a  fully  classical 
description  is  preserved  throughout,  and  leads  to  significant  variances  from  the 
form  of  quantum  calculations. 
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Often  the  initiel  wave functions  and  related  fields  known  as  conjugate 
momenta  (denoted  by  ifr  and  n  in  Figure  1)  are  transformed  to  new  field  variables 
which  are  said  to  diagnonalize  the  Hamiltonian.  In  quantum  systems  these  new 
variables  become  creation  and  annihilation  operators.  In  the  present 
formulation,  they  are  simply  Fourier  coefficients,  and  are  represented  by  n  and 
C  in  Figure  1.  Poisson  bracket  calculations  can  be  made  in  terms  of  the 
diagonalized  fields  as  well  as  the  initial  fields.  In  many  cases  calculations 
are  easier  to  perform  in  terms  of  the  new  diagonalized  field  variables. 

Using  the  Poisson  bracket  relations,  an  equation  of  motion  for  any  field 
variable  can  be  derived  from  the  Hamiltonian.  Many  calculations  are  based  on 
these  Hamiltonian  equations  of  motion,  which  are  fully  equivalent  to  the  origin 
partial  differential  equation  (e.g.,  Navier-Stokes  equation)  to  be  solved.  An 
alternative,  an  extremely  powerful  procedure,  is  to  find  the  Green's  function,  G, 
corresponding  to  H.  If  the  initiel  partial  differential  equation  is  thought  of 
as  an  operator,  0,  operating  on  the  field  variables,  #,  then  G  plays  the  role  of 
the  Inverse  operator  to  0.  The  Green's  function  G  to  the  full  Hamiltonian,  H, 
can  be  derived  from  a  perturbation  expansion  about  the  Green's  function  Gq  to 
the  linearized  problem.  Gq  is  easily  derived  from  the  linearized  Hamiltonian, 

Hq,  so  that  a  procedure  for  solving  a  complex  problem  can  be  built  up 
systematically  beginning  with  the  solution  to  a  simpler  problem.  Green's 
functions  are  central  quantities  in  field  theory.  From  them  the  densities  of 
states,  densities  of  excitations,  transition  probabilities,  lifetimes  of 
excitations,  dispersion  relations,  total  energy,  equation  of  state,  and 
partition  function  can  be  derived. 


Benefits 


One  of  the  hallmarks  of  the  field  theory  approach  is  that  the  development  of 
the  entire  theory  is  completely  systematic,  once  the  true  starting  point,  the 
development  of  a  functional  L  known  as  the  Lagrangian  (from  which  H  can  easily 
be  derived)  has  been  attained.  Thus  theoreticians  skilled  in  field  theory  will 
know  immediately  how  to  perform  relevant  calculations  similar  to  those  done  in 
quantum  systems,  given  the  appropriate  forms  for  L  and  H.  Novices  need  not 
learn  quantum  field  theory  to  perform  calculations.  They  can  begin  at  once  with 
the  rules  for  drawing  and  summing  diagrams.  Indeed,  there  is  a  modern  trend 
toward  teaching  quantum  field  theory  beginning  with  diagrams,  without  any 
reference  to  a  Lagrangian,  or  a  deep  discussion  of  the  theory.  Hence  many-body 
calculations  can  be  made  accessible  to  a  wide  community.  This  is  in  keeping 
with  the  goal  of  the  research  described  in  this  report:  to  facilitate  accurate 
calculations  of  the  properties  of  real  turbulent  flow  fields  in  practical 
systems . 

Notwithstanding  the  main  goal  of  a  practical  calculations!  tool,  it  is  worth 
noting  that  the  Lagrangian  can  be  analyzed  to  obtain  a  deeper  insight  into  both 
the  properties  of  a  many-body  system  and  the  character  of  the  theory.  Given  the 
Lagrangian,  the  stress-energy  tensor,  W  (the  stress-energy  tensor  is  usually 
denoted  by  T  in  quantum  field  theory),  can  be  derived;  and  the  symmetry 
properties  of  L  and  W  can  be  analyzed  to  obtain  all  the  conservation  laws  for 
the  system.  Examples  in  fluid  mechanics  are  the  conservation  of  mass,  energy, 
and  momentum.  Note  that  these  conservation  laws  are  not  postulated  in  addition 
to  the  Navier-Stokes  equation,  but  can  be  systematically  derived  directly  from 
L,  i.e.,  a  collection  of  initial  postulates  is  not  needed  to  describe  the 


system,  only  the  -initial  Lagrangian  and  a  systematic  procedure  for  developing 
the  theory.  This  feature  gives  the  theory  a  strong  foundation. 

Equations  of  notion  for  the  fields  can  be  derived  froa  the  Lagrangian  by 
postulating  that  the  variation  in  the  functional  L  Mill  be  a  minimum,  and  that 
the  fields  are  fixed  at  the  boundaries.  Then  integration  by  parts  and  the 
calculus  of  variations  lead  to  Euler*Lagrange  equations  of  aotion  (see  Appendix 
A).  Hence  the  true  starting  point  for  the  development  of  a  field  theory  should 
be  the  correct  formulation  of  the  Lagrangian.  This  task  has  been  accomplished 
in  the  program  and  is  the  most  important  result  of  the  effort. 

The  value  of  deriving  a  Lagrangian  for  fluid  aotion  has  long  been 
recognized,  and  there  have  been  several  previous  attempts  at  this  task.  Those 
efforts  were  largely  unsuccessful,  primarily  because  of  the  difficulty  involved 
in  accounting  for  dissipation,  but  also  because  of  the  may  in  which  the  pressure 
gradient  term  was  treated.  The  present  theory,  which  we  now  summarize,  departs 
from  earlier  efforts  in  both  respects. 

Summary 

The  technical  description  is  divided  into  three  parts,  which  follow.  In 
Part  I,  the  Lagrangian  is  presented  and  an  analysis  is  made  of  the  quadratic 
terms,  which  correspond  to  the  linear  (noninteracting)  system.  The  Lagrangian 
is  expressed  in  terms  of  the  scalar  and  solenoidal  velocity  potentials,  Y  and 
,  rather  than  directly  in  terms  of  the  velocity.  This  approach  highlights  the 
natural  symmetries  of  the  system  and  brings  in  the  pressure  gradient  term 
without  introducing  an  extra  field.  Velocity  is  accounted  for  by  introducing 
adjoint  fields  for  the  velocity  potential.  These  adjoint  fields  would  grow 
rather  than-  damp  with  time  if  they  were  complex  conjugates  of  the  initial 
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fields.  To  avoid  this  difficulty,  the  adjoint  fields  may  be  considered  to  be 
hypercoaplex. 

It  should  be  noted  that  the  innovation  of  using  hypercomplex  adjoint  fields 
to  describe  terms  of  odd  order  in  the  tiee  and/or  spatial  derivatives  of  an 
equation  of  eotion  is  applicable  to  the  solution  of  a  wide  class  of  partial 
differential  equations,  including  equations  for  passive  scalars  and  many  other 
dissipative  systems. 

Equations  of  motion  for  the  velocity  potential  are  derived  using  Euler- 
Lagrange  equations,  and  found  to  correspond  to  established  results  for  the 
linearized  problem.  The  continuum  mechanics  formalism  is  extended  to  derive  the 
Hamiltonian  and  Poisson  bracket  relations,  as  well  as  the  stress-energy  tensor. 
At  each  step  equations  of  motion  are  generated  anew  and  compared  to  the  initial 
equations  to  check  for  internal  consistency.  The  Hamiltonian  is  then  Fourier 
transformed  and  diagonalized.  Finally,  Green's  functions  for  the  system  are 
derived. 

Part  II  explores  some  of  the  consequences  of  including  the  nonlinear  terms 
(which  correspond  to  the  advective  term  and  higher-order  parts  of  the  pressure 
gradient  term  in  the  Navier-Stokes  equation).  One  advantage  of  the  development 
is  that  the  perturbation  terms  scale  as  the  Mach  number,  rather  than  as  the 
Reynolds  number,  as  is  the  case  in  earlier  works.  However,  since  the  advective 
term  is  of  second  order  in  the  velocity,  the  corresponding  terms  in  the 
Lagrangian  are  of  third  order  in  the  fields.  This  immediately  causes  a  break¬ 
down  of  phase  invariance,  a  problem  which  does  not  arise  in  quantum  field 
theory.  This  may  mean  that  the  Lagrangian  must  be  transformed  to  new  variables. 


The  discussion  in  Part  II  is  basad  on  terms  of  the  Fourier-transformed 
fields.  These  fields  are  further  transformed  to  make  contact  with  various 
calculations  made  on  such  critical  phenomena  as  superfluidity  and 
superconductivity.  Modified  equations  of  motion  are  also  derived  for  the 
original  Fourier-transformed  variables,  and  comparisons  are  made  mith  the 
equations  for  various  critical  systems. 

In  Part  III  procedures  for  making  diagram  calculations  are  indicated  and  a 
simple  illustration  of  a  diagram  calculation  is  given.  An  extensive 
renormalization  calculation  of  the  effective  viscosity,  due  to  John  Erdei,  is 
also  presented.  Next  a  procedure  for  deriving  lagrangians  corresponding  to 
equations' of  motion  for  passive  scalars  is  discussed,  with  an  extensive  example 
of  the  heat  budget  equation.  Justifications  for  several  lengthy  procedural 
steps  are  given  in  the  Appendices. 

Parts  I,  II,  and  III  are  each  written  as  self-contained  works,  complete  with 
their  own  references  for  convenience. 


1 .  INTRODUCTION 


Flowfields  of  practical  importance  are  usually  hard-driven,  nonlinear 
systems  of  great  complexity  that  are  far  from  their  final  relaxed  state  which 
we  term  equilibrium.  Traditional  approaches  to  describing  such  systems  yield 
diminished  results  and  are  not  amenable  to  application  to  a  variety  of  systems 
possessing  diverse  boundary  conditions.  By  contrast,  Many-Body  Theory  is  a 
description  of  a  complex  system  which  has  proved  quite  useful  in  a  variety  of 
areas  that  include  strongly  fluctuating  quantities. 

Experimental  aspects  of  turbulent  flow,  such  as  the  slow  decay  of 
large-scale  structures  under  appropriate  operating  conditions,1  the  bifurcation 

and  quasiperiodic  behavior  of  characteristic  modes  preceding  turbulent 
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behavior,  and  the  onset  of  turbulence  at  critical  Reynolds  numbers,  are 

reminiscent  of  critical  phenomena.  A  comparison  between  turbulence  and  critical 

systems  would  be  greatly  augmented  by  the  use  of  a  Many-Body  formalism.  It 

would  be  especially  useful  for  semi quantitative  calculations  of  scaling 

behavior. 

A  Lagrangian/Hamil ton i an  formalism  has  been  developed  for  a  single¬ 
component,  viscous,  subsonic  flow  in  three  dimensions.  These  functions  will 
describe  the  full  Navier-Stokes  equation  without  approximation.  Results  include 
the  eigenfunctions  of  the  linearized  system  implicitly.  Only  the  fundamental 
bracket  relations  are  required.  Poisson  bracket  equations  are  independent  of 
coordinate  systems.4  Results  can  be  applied  to  various  boundary  conditions  by 
explicitly  inserting  the  particular  eigenfunctions  for  those  conditions. 


Since  the  Hamiltonian  and  the  Poisson  brackets  can  be  shown  to  be 
canonically  invariant,  the  initial  basis  functions  may  readily  be  transformed 
to  other  variables  more  suitable  to  calculation.  Finally,  an  approach  based  on 
Many-Body  techniques  may  make  contact  with  the  vast  literature  on  field  theory. 

In  this  way,  many  useful  results  may  be  obtained  by  a  relatively  small  effort. 

Much  of  the  early  work  on  field  theoretical  calculations  concerning  fluid 

mechanics  was  prior  to  a  maturation  in  field  theory  and  is  expressed  in  a 

sophisticated  format  which  cannot  readily  be  disseminated  to  a  broad  community. 

Attempts  at  applying  diagrammatic  techniques  borrowed  from  quantum  field  theory 

5  6  7 

to  turbulence  in  incompressible  flow  have  been  made  by  Wyld,  Lee,  Edwards, 

go  ^0  i  *i 

Kawasaki,  '  Martin,  et  al.  and  Gledzer  and  Monin.  Except  for  Kawasaki,  who 
used  the  Lagrangian  equation  as  a  model,  all  of  these  authors  simply  made 
perturbation  expansions  directly  from  the  Navier-Stokes  equations,  including  the 
entire  pressure  gradient  term  with  the  advective  term.  Mo  systematic 
development  from  Poisson  brackets  was  considered;  rather  intuition  was  used  to 
renormalize  perturbation  series  that  began  with  the  Reynolds  number  as  an 
expansion  parameter.  It  is  not  surprising  that  these  expansions  tended  to 
become  intractable  after  the  inclusion  of  the  fourth-order  term.  While 
demonstrating  that  the  application  of  field  theory  to  fluid  mechanics  is  a 
plausible  calculation,  and  providing  useful  bases  for  comparison,  this  work  did 
not  follow  a  systematic  development  based  on  a  canonical  formulation  of  the 
governing  partial  differential  (Navier-Stokes)  equation,  as  is  standard  in  field 
theory . 

12  13 

A  Lagrangian  density  corresponding  to  the  Navier-Stokes  equation  '  has 
been  described  in  earlier  work.  However,  this  Lagrangian  (expressed  in  terms  of  . 
the  velocity  itself,  rather  than  its  potential)  is  developed  from  a  Lagrangian 


t  I'l.i'a  t«i  .'^1  «.*  g  * k  A-  M  g  *';  1  *'  .  fl  »<t  g't  .  »  < 


corresponding  to  sound  waves  in  a  Lagrangian  coordinate  system,  augmented  by 
continuity  equations  as  constraints.  Viscous  damping  is  brought  in  by  "posits" 
rather  than  derived  directly  from  a  term  in  the  Lagrangian. 

There  are  many  ways  of  constructing  the  Lagrangian  density  upon  which  the 
Poisson  bracket  relations  are  based.  Several  different  Lagrangians  for  Eulerian 
flow  have  been  described  by  various  authors.14  Host  of  these  are  expressed 
directly  in  terms  of  the  velocity  and  pressure.  The  resulting  formalisms  are 
awkward.  Moreover,  they  avoid  the  truculent  dissipative  term.  The  Lagrangian 
density  described  below  is  equivalent  to  the  full  Navier-Stokes  equation  for 
subsonic  flow  and  is  expressed  in  terms  of  the  scalar  potential  function  and 
the  solenoidal  stream  function.  This  is  a  more  natural  approach,  as  it 
implicitly  couples  the  velocity  and  pressure  and  nicely  separates  potential  flow 
from  rotational  flow.  In  addition,  the  symmetries  and  transformation  properties 
of  the  system  are  more  apparent  using  this  formalism. 


Me  begin  our  formalism  with  a  fluid  density  p,  a  scalar  potential  T,  and  a 
solenoidal  potential , .  Since  the  dissipative  term  in  the  Navier-Stokes 
equation  is  of  even  order  in  the  space  and  time  derivatives,  while  the  other 
terms  are  of  odd  order,  it  is  necessary  to  invoke  adjoint  fields.  This  is  done 
by  expressing  the  Lagrangian  density  in  terms  of  products  of  new  velocity  poten¬ 
tials  and  their  adjoint  fields,  and  then  expressing  T  and  /J  as  the  sum  of  these 
fields.  This  procedure  requires  that  the  "scalar"  velocity  potentials  and  their 
adjoint  fields  become  pseudo-scalars  in  four  spaces  and  that  the  "vector" 
potentials  and  their  adjoint  fields  become  axial  vectors.  It  would  be  nota- 
tionally  simpler  to  work  with  the  complex  conjugates  of  the  potential  fields 
rather  than  with  the  true  adjoint  fields.  This  procedure  would  require  that  we 
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express  odd-order  derivatives  of  Y  and  ^  in  teras  of  the  differences  of  the 
corresponding  derivatives  of  the  new  potential  fields  and  their  adjoint  fields, 
rather  than  their  sums. 


Expressions  for  new  variables  are  obtained  on  the  basis  of  the  syaeetry 
properties  of  the  Lagrangian.  The  Poisson  bracket  relations  for  the  syeeetry- 
generating  variables  are  also  exaained.  It  is  not  necessary  to  use  a  Matrix 
forealise  to  obtain  the  diagonal ization.  As  a  corollary,  expressions  for  the 
conserved  currents  of  the  systea  are  obtained,  and  the  key  continuity  equations 
are  again  generated.  This  work  solidifies  the  basis  for  the  forealise. 


2.  NAVIER-STOKES  LAGRANGIAN 

Me  wish  to  solve  the  Navier-Stokes  equation  for  an  isothereal,  single- 
coeponent  fluid  with  constant  kineeatic  viscosity,  i»: 


v  ♦  v-YV  ♦  vYxYxv  -  DY*v  ♦  HTP  •  F/p  ♦  T/p  , 

P 


where  P  is  the  pressure,  p  is  the  fluid  density,  and 


0  ■  C/p  +  jjV  . 


(2.1) 


(2.2) 


To  do  so,  we  construct  the  appropriate  Lagrangian,  L,  in  teras  of  fields  for  the 
velocity  potential.  L  is  an  integral  over  a  Lagrangian  density 


L  »  /dtd»r£ 


(2.3) 


Me  express  v  in  teres  of  a  scalar  potential  field  Y  and  a  solenoid  potential 


field  A 


v  »  YY  ♦  Vxrf  , 


(2.4) 


l 


V.V.v.v.v.v;.- vv:/: 


with 


A  »  A  ♦  A  (2.5) 

and 

*«♦♦♦,  (2.6) 

where  ~  denotes  the  hypercoaplex  adjoint.  We  sake  the  approximation 

P  -  -  pf  -  \p  [  lVT|2  -  f2  ]  ♦  PQ  ,  (2.7) 

c 

where  PQ  -js  a  constant.  ‘  Higher-order  approximations  for  P  may  be  obtained  by 
adding  correction  terms  to  Equation  (2.7).  Denoting  the  speed  of  sound  by  c, 
the  required  Lagrahgian  density  is  given  by 

£(+.  *  ,  A,  A  )  *  £^  *  £k  *  £t+  £hl  ♦  £p  ,  (2.8) 

where 

2 

2“j^  ■  it  -  ^({YqV2*  -  ~  c2**.**  +  ♦  p)£,  (2.9) 

2  # 

2^2a  -  A»A  +  ^v(Ayo»VxVxA  -  7x7xAyo*A)  ♦  1(A  +  A)B,  (2.10) 


®p^I  “  2  +  *V<VxA  "  7x^),(V^  ”  **) 

♦  2(Yo*  *  *V<7xA  "  7x>i),(7xA  -  VxA) , 
c2  .1 

tv  (r0*  ♦  *r0)  <*►  -  **)•(**  -  tip)  , 

2 

2p  4  •  <v +  [**  -  ^i2  ]  -  V*  -  i)2 . 

c 


(2.11) 

(2.12) 

(2.13) 
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In  Equation  (2.9),  wo  let 


F  »  -  VU 


(2.14) 


Similarly,  we  let 


e  ■  -  U/p  . 


B  -  (Vxt)/p  . 


(2.15) 


(2.16) 


£^  describes  the  potential  velocity  field,  £^  describes  the  solenoidal 
velocity  field,  and  £^  describes  the  interaction  between  these  two  fields  due  to 
the  advective  tera  in  the  Navier-Stokes  equation.  The  nonlinearity  in  the 
Navier-Stokes  equation  is  contained  in  £^  and  ^NL  and  in  correction  terms  to 
Equation  (2.7). 


2.1  Equations  of  Motion 

The  Euler- Lagrange  equations  of  motion  for  ♦,  A,  and  A  are  obtained  by 
minimizing  L  with  respect  to  these  fields  through  the  calculus  of  variations, 
using  integration  by  parts  (see  Appendix  A).  These  equations  may  then  be  used 
in  conjunction  with  Equations  (2.4)  through  (2.6)  to  obtain  an  equation  of 
motion  for  v.  Since  the  potential  field  equations  will  be  of  second  order  in 
time,  the  resulting  equation  must  be  integrated  with  respect  to  time  to  regain 
Equation  (2.1).  Calculations  based  on  ^will  be  made  by  treeting  £  ,  £^  and 
the  nonlinear  terms  in  the  series  approximation  for  P  as  perturbations  about  the 
lineer  tera  in  £^  and  £^.  In  general  the  expansion  parameters  should  go  roughly 
as  *(v/c),  so  the  series  should  converge  rapidly  for  subsonic  flow.  As  the 
development  of  the  formalism  proceeds,  the  formelism  will  take  on  a  structure 
that  closely  parallels  that  of  quantum  field  theory.  Nevertheless,  the  system 


& 


i.rwMJ 


remains  completely  classical,  as  explained  in  Appendix  B.  Me  begin  by 
investigating  the  linear  Lagrangian  density^,  where 

£Q.£*+£^  *  (2‘17) 

Variation  of  9  in  £Q  gives 

••  2*  2  2  1 

♦  -  DrovS  -  cT*  -  ~e  .  (2.18) 

A  corresponding  variation  of  in  £Q  gives 

et  e 

*  ♦  DV2^Y  -  C2V2$  »  -£  .  (2.19) 

o  z 

Variation  of  A  in  £  yields 

A  ♦  VyJxVxk  ■  ,  (2.20) 

while  variation  of  A  In  £Q  results  in 

-  It 

A  -  vVxVxAyq  -  2s  (2.21) 

These  results  are  in  agreement  with  those  of  other  authors. C'5'16)  -me  momenta 
conjugate  to  p,  <Jr,  A,  and  A  are  respectively: 

ir  ■  a  £  /a*  -  -^(i  *j0v2*rj.  '(2.22) 

0  2C2  2  0 

s.ai0/aj-  ^2($  -  50rov2#).  (2.23) 

p  ■  a£0/aA  »  -^(A  -  ^vfxVxAr0) , 


(2.24) 


and 


p  •  a  -£q/3A  ■  -%(A  ♦  |wy07x7xA). 
2c  v 


2.2  Haniltonian 

The  Haniltonian  density  corresponding  to  Jt ^  is  given  by 


-2-14 

2c2 


[i*  *  c2vm>  - 1(*  +  ♦)«] 


2c 


r  2c2  1  2- 

[_("7  "  '  2  DV  ^ 


0)(f^  *  Kv^> +  c 


-  +  ^>e 


■] 


Sieilarily,  the  Hamiltonian  density^,  corresponding  to^  is 


given  by 


P*A  +  p*A  -  £a 


2c 


j^A*A  -  |(A  ♦  A)bJ 


2c 


P 


+  P  -  KVX7XA) 


-  |(A  ♦  A)bJ 


The  Haniltonian  equations  of  notion  for  n,  < (r,  and  n  are 


<fr  -  «H/3w  »  ^  ir  ♦  ^Oy/2*, 


n  -  -  -er(-  jOT2#*  ♦  c272*  +56), 

2c  ‘ 
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(2.25) 


(2.26) 


(2.27) 


(2.28) 


(2.29) 


(2.30) 


(2.31) 


a 


a 
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Li 


$ 
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•Sb 


M 


ai 


7x 4-  y^VxA  -  VxAy^.  ( 

Then,  adding  Equations  (2.18)  and  (2.19),  we  obtain 

..  2*  2  2 

7  -  07  7  -  c  7  t  =  e  ( 

Adding  Equation  (2.20)  and  (2.21),  we  find 

/}  +  v7x7x/?  =  B  ( 

2.3  Nonlinear  Terws 
Assuming  the  gauge  equation 

7*A  =  0  ,  ( 

the  addition  of  the  second  expression  in  to  does  not  affect  Equations 
(2.20)  and  (2.21),  but  adds  the  term 

iTof;[<7xA  '  ?xA)-(VxA  -  7xA) ] 

to  both  Equation  (2. in)  and  Equation  (2.19).  The  corresponding  change  in 
Equation  (2.40)  is 

|^(7xrf-7xrf)  . 

The  addition  of  L...  to  L.  adds  the  terms 
NL  i p 

-{yjf  *  bQ)  (*2*  -  72*) 

on  the  left-hand  side  of  Equation  (2.18)  and 

(Ttp  -  7»j»)  •  (y  Tijf  ♦  )  ♦  o  (T  'l'  ♦  'bj(V2'\>  ~  72tj>) 

0  0  c  0  0 


[20) 


.  *•"«  V* .  "*- *■"*-  •"  _  «*« 


i] 


on  the  left-hand  side  of  Equation  (2.19).  The  corresponding  change  in  Equation 
(2.42)  is 


(V*  -  VjO*(Y0W>  -  tJy0) 

«  |d(VT«V¥)/at 

Me  may  integrate  Equations  (2.42)  and  (2.43)  with  respect  to  time,  and  use 
Equations  (2.4),  (2.15),  and  (2.16)  to  obtain  the  linearized  counterpart  to 
Equation  (2.1). 

The  addition  of  the  first  expression  in  £^  to  £Q  does  not  affect. Equations 
(2.20)  and  (2.21),  but  adds  the  terms 


Yolt  [(7*  ‘  W*<VxA  ~  V*A) 


+  (Y0*|»  +  Vihro)»(VxA  -  VxA) 


to  the  left  side  of  Equation  (2.18)  and  the  terms 


To!*  [(**  -  **).(VxA  ~  VxA)]  -  (yoV*  +  v4yo).(VxA  -  VxA) 

to  the  left  side  of  Equation  (2.19).  The  corresponding  additions  to  Equation 
(2.42)  are  terms  of  the  form 

(Vx^j  +  ^W.(Vx^), 

provided  again  that  Equations  (2.36)  and  (2.37)  hold.  As  a  check,  we  differ¬ 
entiate  Equation  (2.1)  with  respect  to  time,  and  substitute  in  for  v  using  the 
Equation  (2.4).  We  make  use  of  the  linear  approximation  for  the  pressure 
gradient  term, 

<P  +  V»v  »  0,  (2.45) 


[21] 


to  obtain  the  relations 


Vx| fl  +  vVxVw<7  ]  -  8  ■  0 


(2.46 


and 


V[¥  -  0 vH  -  l-p2*  ♦  |  |^(TT.?¥)  ♦  |^(7>W  *Tf )  +  |^(Vx^‘?xtf)  -  e]  •  0.  (2.47) 
c 


2.4  Poisson  Bracket  Relations 

Me  will  define  Poisson  brackets  in  teres  of  a  variational  derivative: 


[X,Y] 


m  r  ax  ay  ax  by 
"  ^  1  an4  "  an{  «*i 


(2.48 


where  the  fields  #.  are  given  by  7,  f  ,  A,  and  A  and  the  n  fields  are  their 

•  i 

conjugate  aoeenta  n,  n,  p,  and  p. 


Me  have  the  usual  relations 

[♦.j (r) ,  * j ( r ' ) ]  -  [n^(r) ,  n^r')]  -  0.  (2.49 

For  the  conjugate  fields  we  write 

C*i  (r) ,  lyr1)]  -  C  ^  fl(r-r').  (2.50 

For  the  scalar  fields  we  have  C  *  1,  but  for  the  vector  fields  we  have  C  ■  a: 
we  shall  show  in  Section  3. 


Keeping  in  mind  that  our  Poisson  brackets  are  defined  in  terms  of  varia¬ 
tional  derivatives,  the  equations  of  notion  expressed  in  terms  of  Poisson 
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brackets  are  in  standard  form: 


[.>  m  LLi 


p  *  [p»  H]  +  ft  ' 


(2.53) 


for  any  variable  F.  Me  can  check  our  formalism  for  consistency  by  using  the 
relation  (2.53)  to  retrieve  our  Equations  of  motion.  With  given  by: 

"a  ■  ^2  iryv  *  J»7x»xAr0  j  •  .  I^VxTka)  (2.54) 


we  find 


A  »  [A,  Ha] 


2 [A.  p]  p  -  ^vyo7x7xA 


2£  p  -  1wt07x»xa  . 


(2.55) 


(2.56) 


Similarly,  recalling  that  our  brackets  are  defined  in  terms  of  variational 


derivatives. 


P  3  TP'V 


i_  +  JvVxVxAyoj  ^VYQ[VxVxA,  p] 


~^~zvy  7x7x  f^-p  +  ^i>7x7xAy 
._2  o  V  P  2  'c 


-firV7x7xA 

4c 


(2.57) 


(2.58) 


/v*/  -/o  • 


With  given  by 


2  - 
c  V* 


•  V* 


( 


we  have 


1  2 
+  ^0yQ7zKf 


( 


Similarly, 


*  *  C»r.H^3 

2 

*  '  2°v2*0  'fa*''72*]  *  c2W-[x,  ViJ»]  ] 

2c 

■  ^  ['  i07'  (¥*  -  K^o)  +  cV*  ] 

*  *^2  (-  ^2ir0  +  cV*) 

2c 


in  the  absence  of  external  forces.  In  the  same  manner  we  can  reproduce  the 
equations  of  motion  for 


3.  CONSERVED  CURRENTS 

The  prediction  of  conserved  currents  on  the  basis  of  the  symmetry  of  the 
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Lagrangian  density  £  is  discussed  in  detail  by  Roman  and  by  Handel  and  Shaw. 
Although  their  discussion  is  directed  towards  the  quantum  systems  of  particle 
physics,  all  the  arguments  can  be  applied  to  a  classical  system.  Nothier's 
theorem  states  that  the  invariance  of  the  Lagrangian  density  with  respect  to  a 
symmetry  transformation  implies  the  existence  of  a  conserved  variable.  Consider 

a  four-vector,  f^,  such  that 

a^f"  *  o.  (3.i) 

Integrating  Equation  (3.1)  over  coordinate  space,  and  making  the  identification 

F  •  /d3rfQ  ,  (3.2) 


yields 


3 


-  /d3rX?Va*f  . 

(3.3) 

-  JdAf  -  0 

(3.4) 

Consequently,  if  the  variation  of  £  with  respect  to  some  symmetry  trans¬ 
formation  yields  an  equation  of  the  form  of  Equation  (3.1)  for  some  four  vector 
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f^,  we  are  led  to  a  conserved  variable  F.  As  shown  in  Handel  and  Shaw,  gauge 
invariance,  expressed  by 


*  -  e16<fr 

(3.5) 

$  -  e~ie*  , 


leads  to  the  conservation  of  the  quantity 

/d  rfx#  -  n  $  ]  =  m/p 


(3.6) 


[25] 


The  continuity  equation  for  m  is 


m  ♦  7»j  ■  0 

with 

j  -  *7*  -  #7$ 

Invariance  with  respect  to  the  four-translation 


<fr(x^)  -  (x^) 

A(xy)  -  A' (x^) 

leads  to  the  conservation  of  the  stress-energy  tensors  and 
(see  Appendix  0). 

The  general  fore  of  the  eleeents  of  T*  is  given  by: 


T* 

»  df  — + 

2lU 

VO 

V**a 

vp*op 

♦  #^  + 

2l[^ 

vd*o 

pi  Vptyap 

(3.7) 


(3.8) 


(3:9) 


(3.1 


This  definition  of  is  a  necessary  extension  of  the  stress  tensors  described 
by  Horse  and  Feshback16  and  by  Morse  and  Ingard15  if  dissipation  is  to  be 
included  in  L.  A  sieilar  definition  will  hold  for  W*.  Explicit  forms  for  the 
eleeents  of  the  stress  tensor  are 


+  c27$7>|»  -  1(*  +  <|>)e 
2c  *•  Z 


H 


(3.11 


lr 


2c 


-lD(7^mroi|>  -  \lnrQ Vty) 


c  (t|f7«*r  -  4f7t{») 


(3.12 
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^  v.  /.  o'.  -  ^  >  ,  »•.  .  .  -  .  *  .  1  .  *  . 


T  *  +  ir9^)  ■  -^rP  (3.13) 

2c£  2c 

Txy--?*x*y  '  •**•  *  (3-,4> 

and 

Txx  *  ^[**  *  “V****  -  **xx>  -  e*<*yV  *z*z  -  W  •  '"=•]•  <3-,5> 

The  components  of  satisfy  the  continuity  relations 

H  +  7*1  •  0  (3.16) 

and  P  ♦  V«W  *  0  ,  (3.17) 

where  W  is  a  3x3  tensor  containing  the  spatial  parts  of  T^,  known  as  the  wave- 
stress  tensor . 

We  see  that  the  continuity  equation  for  the  stress-energy  tensor  follows 
directly  froe  the  fore  of  our  Lagrangian  density,  i.e.,  we  can  derive  it  from 
our  Lagrangian;  it  does  not  represent  a  subsidiary  condition  on  our  formalise. 

Although  £q  is  invariant  with  respect  to  gauge  transformations,  the 
nonlinear  terms  in  our  full  Lagrangian  density  break  gauge  symmetry.  Neither 
the  A  nor  the  *  fields  can  be  aade  to  satisfy  local  gauge  invariance  in  any 
simple  or  reasonable  manner.  However,  the  full  Lagrangian  density  is  preserved 
with  respect  to  four- translations. 

The  symmetries  of  the  Lagrangian/Hamlltonian  density  are  important  not  only 
for  generating  additional  equations  of  motion,  but  because  they  provide  valuable 
clues  and,  in  fact,  a  systematic  procedure  for  examining  the  behavior  of  the 
system  when  it  undergoes  a  change  of  state,  e.g.,  from  laminar  to  turbulent 

#1  nu 


Since  the  conserved  variables  may  be  derived  from  the  symmetry  properties  of  the 
Lagrangian,  their  Poisson  bracket  relations  with  other  variables  offer  clues  to 
the  behavior  of  the  system  with  respect  to  various  symmetries.  In  particular, 
if  F  is  a  conserved  variable,  and  if  there  exist  (possibly  the  same)  variables 
such  that 

[F(0),  X]  *  Y  (3.18) 

Mith 

<Y>  *  0,Y  $  0  (3.19) 

c 

for  some  parameter  8  and  a  critical  value  of  that  parameter,  0  then  Equations 
(3.18)  and  (3.19)  signal  a  breaking  of  symmetry,  usually  associated  with  a  change 
in  state,  at  8  .  Me  seek  such  relations  in  an  attempt  to  develop  a  theory  of 

G 

coherent  structures  at  the  onset  of  turbulence.  Consequently,  it  is  of  interest 
to  examine  various  Poisson  bracket  relations  pertaining  to  conserved  variables. 

Me  have: 


lm.  *] 

«  -  (#  , 

(3.20) 

[■.  *] 

•  pi 

(3.21) 

[a,  x] 

•  pn 

(3.22) 

[»,  x] 

■  -  px 

(3.23) 
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Equations  (3.1)  and  (3.10)  to  (3.13)  lead  to  the  conservation  of  P  in  the  form  ' 

P  ■  /d3r(x9ifr  ♦  xV*]  (3.24) 

The  function  P  has  the  property  that: 

[F,  P]  ■  7F  .  (3.25) 

From  Equations  (3.18)  and  (3.7)  we  see  thet  if: 

<v>  ■  <V*  -  4  0  (3.26) 

then  translation  symmetry,  generated  by  P,  will  be  broken. 


4. 


DIAGONAL I ZATION 


We  Mill  now  make  linear  transformations  of  the  scalar  potentials  which  will 
diagonalize  in  the  absence  of  external  driving  forces.  It  will  be  helpful  in 
performing  the  required  algebra  to  define  a  frequency  given  by 


“k 


c2k2  +  Jo2k4, 


(4.1) 


for  a  mode  with  wavenumber  k.  We  will  ignore  external  driving  forces  in  the 
following  discussion. 


We  first  expand  +  and  ♦  in  terms  of  the  normal  modes  of  system.  We 
express  the  normal  modes  as  functions  of  their  wavevectors  k  through  the 
argument  ik*r.  We  will  assume  a  stationary  system,  and  express  ^(t)  in  terms 
of  the  phases 


■  e 


'v* 


(*•2) 


■  e 


t 


(4.3) 


*♦  . ,  V 

8Ck  *  * 


and 


■  e 


t 


(4.4) 


(4.5) 


The  scalar  fields  are  expanded  as 


where  we  use  the  notation 


0k  -  u*yA  (4.8) 

The  constant  factor  c/Vpu^  is  inserted  into  the  definitions  given  by 

Equations  (4.6)  and  (4.7)  to  produce  expressions  for  quantities  of  interest  in 

teres  of  the  Fourier  transforaed  fields  that  are  compact  and  that  highlight  the 
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analogy  with  work  in  other  areas  using  second  quantization.  The  system  is 

"quantized"  only  to  the  extent  that  we  are  using  expansions  in  Equations  (4.6) 
and  (4.7)  in  series  rather  than  continuous  transforms.  Our  fields  are  classical 
fields,  and  we  have  not  as  yet  introduced  eigenvectors  or  operators. 

Two  species  of  functions,  nk  and  (k  in  the  transformations  given  by 

Equations  (4.6)  and  (4.7)  are  needed  because  the  equation  of  motion  for  t|r 

contains  terms  that  are  of  zero,  first,  and  second  order  in  time,  as  well  as 

terms  that  are  of  third  order  in  the  field  strength.  The  presence  of  terms  of 

odd  order  in  either  the  time  or  space  derivatives  or  in  the  field  strength 

dictates  the  use  of  complex  fields.  If  the  equation  of  motion  is  genuinely  of 

second  order  with  respect  to  time  and  the  Lagrangian  is  expressed  in  terms  of 

20 

complex  fields,  then  two  species  are  required  for  diagonal ization.  This 

18—20  24 

occurs  for  example,  in  the  description  of  a  charged  meson  field.  '  If  the 

equation  of  motion  is  of  second  order  and  only  a  real  field  is  required  in  the 

Lagrangian,  then  only  one  species  of  function  is  required  for  diagonal ization. 

19  24—26 

Examples  are  phonon  fields  and  the  neutral  ir  meson  fields.  '  A  second- 

order  equation  is  equivalent  to  two  first-order  equations,  so  complex  fields 
with  first-order  equations  of  motion  can  be  expressed  in  terms  of  a  real  field 


20 

obeying  a  second-order  equation  of  motion..  This  occurs,  for  example,  with 
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the  Schrodinger  equation  '  and  the  diffusion  equation.  The  equation 
of  motion  for  the  A  fields  is  similar  in  form  to  the  latter  two  cases,  and  we 
will  need  only  one  species  of  function  to  diagonalize  H^,  as  we  shall  show 
below. 

We  emphasize  that  Equations  (4.6)  to  (4.8)  show  how  the  development  may  be 
applied  to  any  geometry.  In  the  remainder  of  the  discussion  we  will  concentrate 
on  the  functions  ^(t)  and  $^(t)  and  their  derivatives.  This  discussion  will 
be  completely  general.  To  apply  the  results  to  a  given  geometry,  we  need  only 
substitute  the  appropriate  functions  u^ik^r)  and  Q^t-ik'r).  In  the  discussion, 

whenever  we  refer  to  these  functions  explicitly,  we  will  use  the  plane  wave 

ik*r  — ik*r  -  — ik*r 

solutions  u^e  and  u^yAe  ■  u^e 

4.1  Scalar  Potential 
Since 

5  -  0yo72*  -  cV*  »  0  ,  (4.9) 

and 

*  +  0V2Jy  -  c2V2iji  »  0  ,  (4.10) 

o 


we  have 


2  2  2  2 
V  *  iDyoV  - c  k 


(4.11) 


and 


2  2  2 

iDy  u,.k  -  c  k  ■  0 
o  Ck 


(4.12) 


Equations  (4.11)  and  (4.12)  imply  that 


Note  that 


“nk  “ck  *  c2k2  ' 


so  that  substituting  Equations  (4.6)  and  (4.7)  into  the  relation 


Jt*  -  -4r  [£i  +  c27*.V*] 


♦  2c2 


yields 


H*  *  I  [(v  * cV)  Vk  *  ("{k  *  oV)<ck]  / 


2  2 
+  c  k 


2\  • 


assuming  that 


/dr  uk(r)  uk,(r)  »  «k  k, 


i.e.,  that  the  ufc  are  orthonormal. 


We  now  rewrite  Equation  (4.11)  in  the  form 


u5k  *  °2k2  ■  2(v  *  w<v2>v  ' 


or 


V  *  °2r2  ■  2“k“,k 


Similarly,  we  may  rewrite  Equation  (4.12)  as 


Substituting  Equations  (4.20)  and  (4.22)  into  Equation  (4.17),  we  obtain 


at  last 


•v  •  l  ('A  *  "ck#k ) 
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in  standard  fora,  '  1  albeit  with  coaplex  eigenfrequencies. 

be  convenient  to  use  the  coapact  notation 

<Mr,t)  »  £  <Pk(t)  uk(ik*p)  • 


*(r, t)  a  Y  $k(t)  uk(-ik»r)  , 


(4.23) 


It  will  often 


where 


*k  *  (^)5<  <£jk +  \V  ' 


♦k  •  <^)*  *  Vnk}  ' 


(4.24) 


(4.25) 


(4.26) 


(4.27) 


The  fields  ir  and  ir  aay  also  be  expanded  in  a  aanner  similar  to  $  and  tp: 


*  *  I  jrk(t)uk(-ik»r) 


From  the  relations 


*  ■  I  *k(t)uk(ikT) 


n  -  (i  ♦  H0V2*v  ) 


n  »  (*  -  *Dr  V2#) 

2c 


(4.28) 


(4.29) 


(4.30) 


(4.31) 


X VO 


and  Equations  (4.6)  and  (4.7)  we  find 


V*> 


[<%<  '  9Jk  -  (i“ck  *  ,,Dl'ok2,Cj;9Ck 


V*»  *  ^<»/“k>’<  [<-%k  *  Wl'ol'2)"keSk  *  (1"{k  *  ^ok2)Cg0Ck. 


Now  from  Equation  (4.14)  note  that 


V  -  ’“V  m  'uk  *  iuck  *  ’“’V 


consequently  Equations  (4.31)  and  (4.32)  may  be  simplified  to 


V'l  ■  '(pnc)’*/2'  t"k<k  -  <secki 


V*»  ■  i«"k>’,/2c  '^9ck  -  W 


Equations  (4.6),  (4.7),  (4.35),  and  (4.36)  are  the  building  block 

expressions  for  the  remaining  development.  We  now  have  all  the  information  we 

*  * 

need  to  proceed.  First  expressions  for  the  quantities  n^,  n^,  Ck»  and  Ck 
(which  are  independent  of  r  and  t)  are  obtained  in  terms  of  the  original 
fields: 


[k'pVc  *k  *  2c/'/(,“k”k]  9nk  ■ 


\  '  1‘[,'5VC  *k  -  2c'^»k  ]  9Jk  • 


c*  ■  \W^/c  +  2c /i/pujr,,  1 


/  ,• 


v*.  _W*. 


^TSSvS  >»*•**■*>  «-4 


•a 


wvtwwwi«OTvwr> 


C  ■  »[l^/o  *k  -  Zc/l/SyrJ  9*k  . 


(4.40) 


Of  course,  we  aay  also  write 


<-i/*[  ?«• 


n  1  VUk  '  (4>41) 


•  f  33r  ^  *  -  ^=;  I  «V 

2y  c  i/pwj^  nk 


Jk  ' 


(4.42) 


■*-*/* 


'pw^-  *  9CkUk 


(4.43) 


ck  ■  i/d3r  [  *  1 9ck= 


(4.44) 


4.2  Vector  Potential 


The  Hamiltonian  for  the  A  fields,  H^,  is  already  diagonalized,  as  can  be 


seen  from  the  equation 


Jf.  *  -®r  A»A 
A  2C2 


(4.45) 


in  the  absence  of  external  forces.  The  A  fields  are  expanded  as  follows: 


A  -  I  Ak(t)uR(r) 


l  vCa 


(4.46) 


A  -  I  Ak(t)uk(r) 


l  *;e.kGk 


iVIVWWp 


(4.47) 


where 


_  iu  t 

8ak  -  e  a  . 


«.<8) 


and 


(4.49) 


From  the  equations  of  Motion  for  A  and  A  in  the  absence  of  external  forces, 


A  ♦  vy  VxVxA  ■  0  . 
o 


and 


(4.50) 


A  -  ulxVxAy  •  0 
o 


(4.51) 


we  obtain  the  dispersion  relation 


'".k  ■  -"V 


or 


(4.52) 


“ak- 


(4.53) 


in  the  agreement  with  the  well-known  results. 
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Substituting  Equations  (4.46),  (4.47)  in  Equation  (4.45),  we  find 


HA  *  l  ".k  *k\ 


(4.54) 


Vie  may  also  expand  the  momenta  conjugate  to  the  A  fields: 

P  *  I  Pk(t)uR(r) 


(4.55) 


P  ■  I  Pk ( t ) p ) 


(4.56) 


Then,  from  the  relations 


p  «  ,  (A  -  HWxVxAy  ) 

2c 


(4.57) 


(4.58) 


P  *  (A  +  JfVyVxVxA ) 

2c2  0 

mw  find,  using  Equations  (4.46),  (4.47),  (4.52),  and  (4.53), 

pk  *  i^P“ak/8cz  a£  0ak  (4.59) 

and 

K  ’  '  i/p“ak/8c"  \  9Jk  '  (4'60' 

The  inverse  transformations  corresponding  to  H^are 

■k \  ®.k  ■  <4-6” 

<  -  \  ®Jk  •  <‘-62» 

•k  "  ,/BcJ/w.k  ^k  ®ak  '  (4-63) 

and 

<  ■  -"'•^/"‘Wk  0*k  .  (4.64, 

We  preserve  the  expression  for  in  terns  of  the  functions  a£ak  in  standard 
fora,  while  maintaining  Poisson  bracket  relations  in  the  diagonal  representation 
sinilar  to  the  ones  conmonly  found  in  field  theory.  We  recall  that  the  complex 
first-order  equations  for  the  A  fields  (these  fields  are  not  explicitly  of  first 
order,  but  clearly  are  equivalent  to  first-order  equations)  nay  be  related  to  a 
real  field  that  obeys  a  truly  second-order  equation  of  motion.  Let  this  field 


be  defined  by  the  equations 


*k  ■  \  *  V’’ 


Km\-  V1", 


ak  ' 


where 


Then  we  have 


\  *  "ak  °k  *  0 


*k  *  ’"ak*k 


*k  ■  -1uak\ 


Using  our  dispersion  relation  for  the  vector  fields. 


“ak  '  ’  Vs  • 


and  the  relation 


VxVxA^  »  -  k4A^ 


(4.65) 


(4.66) 


(4.67) 


(4.68) 


(4.69) 


(4.70) 


(4.71) 


we  see  that  Equations  (4.68)  and  (4.69)  are  equivalent  to  our  equations  of 
Motion  for  the  vector  fields.  At  the  same  time,  Equation  (4.67)  is  the  equation 
of  motion  for  a  harmonic  oscillator  with  coordinate  o^,  frequency  u^,  and 
conjugate  momentum  given  by 


"ok  ■«k 


(4.72) 


Consequently,  we  must  have 


[V  V3  3  CV  *V3  3  ak,k’ 


(4.73) 


[38] 


tW1  “  <V  V3  a  0 


(4.74) 


Equations  (4.75)  together  with  Equations  (4.71)  to  (4.74)  imply  that 

[V  pk'3  *  CV  V3  "  2akfk'  '  (4.75) 

and 

v>  ■  i4k,k‘  •  (4-76’ 

as  asserted  following  Equations  (2.47). 

Froa  Equations  (3.48)  to  (3.51),  (3.68)  to  (3.71),  (2.45),  and  (2.46)  we 
obtain  the  relations 

[XV  Xj ]  -  [X*.  Xj]  -  o  •  (4.77) 

and 

[X*.  Xj]  -  .  (4.78) 

where  the  fields  x..  are  given  by  )?k  and  f^.  We  aay  show  this  explicitly  as 
follows: 

tnk'  nk'3  s  4(2i[^k'  V1  "  21 CV  *k3)  *  ifik,k*  (4,79) 

and 

[{J,  Ck,3  *  J(2i[*k,  *k.]  -  2iC\'  $k,l)  *  i6k,k'  •  (4*80) 

The  factor  i  in  Equation  (4.78)  is  due  to  our  use  of  Poisson  brackets  to  describe 
classical  fields  rather  than  coaaautators  to  describe  quantua  fields. 


4.3  Conserved  Currents 

Me  can  define  variables  analogous  to  the  number  operators  in  quantum 
mechanics.  Because  we  have  a  classic  system  described  by  Poisson  brackets,  we 
must  introduce  a  factor  of  i  in  our  definitions; 


Nn  ‘  l  V  ■  £  'VV  • 

(4.81) 

"c  *  £  \k  *  ■ 

(4.82) 

and 

"*  ’  k  '  k  ' 

(4.83) 

Mith 

these  definitions  we  obtain  the  standard  relations 

trV'  V3  *  “  nkakk* 

(4.84) 

tnrjk'  nk'3  *  nk*kk' 

(4.85) 

[nCk'  Ck’3  *  ‘  Ckdkk’ 

(4.86) 

[nCk'  fk'3  "  Ck5kk' 

(4.87) 

fnak'  *k'3  Vkk' 

(4.88) 

and 

tn.k-  ’ll'1  *  *kakk’  . 

(4.89) 

Hq  »  W  *  (-*-90) 


where 


and 


We  now  look  for  variables  that  are  conserved  with  respect  to  Hq.  some 
immediate  examples  are  N^,  N^,  N#  and 

*  *  ip^(nk\  '  "  <4.96) 

and 

P  "  ^k(nk\  *  CkCk)  (4,97) 

5.  GREEN'S  FUNCTIONS 

The  Green's  function  Gu(r|r')  is  the  solution  to  the  equation 

<Na  -  w)Qw(r|r')  -  -  «(r  -  r')  .  (5.1) 

In  Equation  (5.1)  we  will  assume  that  HA  is  a  scalar,  as  we  have  previously. 

All  of  the  basis  Hamiltonians  which  we  have  discussed  are  given  in  terms  of  a 
quadratic  product  of  the  field  amplitudes,  which  are  vectors.  Thus  far,  we  have 


[41] 


assumed  that  this  product  is  a  dot  product,  i.e.. 


Hj  ~  <X  |X> 


(5.2) 


for  any  Hamiltonian  H ^  in  teras  of  tha  appropriate  fields  x  and  x 


5.1  Derivations 

6u(r|r')  is  first  expanded  in  teras  of  the  eigenaodes  of  H^,  as  is  coaaonly 
done.  Usually  a  linear  expansion  is  aade,  but  in  anticipation  that 

since  G  is  the  inverse  of  H.  it  will  be  bilinear  in  the  fields,  we  write 

U  A 


„  ,  i  , ,  rr-  _  *  -ik’^r*.  ik»r 

Vr'p  >  -g.W* 


V 


(5.3) 


Substituting  Equation  (5.3)  into  Equation  (5.1)  we  find 
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..  rr  «  /  v  *  -ik'»r'  ik»r 

«(r  -  r»)  ■  -  wak  *k'*  *ke 


(5.4) 


Me  will  continue  with  this  foraulation.  Me  could  also  express  the 
Haailtonians  in  teras  of  tha  vector  product,  |x><X  |,  (See  Appendix  C),  in  which 
case  we  would  have  an  operator  formulation.  In  that  case  the  discussion  would 
parallel  the  dot  product  case,  and  all  of  the  following  treatment  would  apply. 
Most  of  the  earlier  discussion  would  also  be  valid,  but  it  would  be  necessary 
to  replace  the  Poisson  brackets  with  coaautators. 


Since 


.  *<-  >!-  tl.  tl.  il. 


i.  ft.  |V*t . 


Jta  |V| 


in  Equation  (5.10)  rather  than  summing  over  k  gives 

G(k,t-t' )  »  G(k,T) 

-  il  a-aV^ 
k  *  * 


(5.15) 


■  dr  G(r,t|r* ,t' )  .  (5.16) 

Proceeding  similarly  Me  may  find  the  Green's  function  S  (r|r')  corresponding  to 


H. .  S  must  satisfy  the  equation. 

f  (a) 


(H^  -  u)Su(r|r')  =  -<5(r  -  r') 


(5.17) 


Assume18,19  that  S  is  of  the  form 

(l) 


Then,  substituting  Equation  (5.18)  into  Equation  (5.17),  we  find,  from  Equation 
(4.23)  and  an  argument  parallel  to  that  above,  that 


SH  .  — SSL. 

akk'  u>  -  u 


(5.19) 


Therefore  Me  may  Mrite 


V|r’> 


u  -  <a. 


«*')  -  S' (r 

w 


(5.20) 


(5.21) 


I  H5-5- 

k  1  «  •  “, 


ik-(r-r' ) 


(5.22) 


.’Xv, .  \v 


The  Fourier  transform  of  S^,  S(r,t|r’ ,t' ) ,  is  then  given  by 
S(r,t|r',t')  *  9(T)Sn(r,tlr' ,t’ )  -  0(-r)Sf (r,t|r* ,t* ) 


Jl  rfjVllS - VK_L,«f(i-' 

J2*  k^“  -  V  “*"4 


r)  -  iw(t-t') 


HlfdJ— ! - 

J  U>  -  “  *  “,  J 


“nk  ’  “k  '  iDl,ok 


V  '  “k  *  il”.1' 


the  contour  integral  gives 


Sn(r,t|r',t')  -  i«(r  -  r' )e“1ukT“ok  T  ,  t  >  f 


0  ,  t  <  t\ 


SC(r,t|r\t')  *  -i«(r  -  r'Je’V-0*  T  ,  t  <  t' 


■  0  ,  t  >  t'  . 


From  Equations  (5.25)  to  (5.27)  and  (4.15)  m  find 


S(r,t|r' ,t' ) 


-iw(t-t') 


wz  -  ♦  2iwrQDkk 


Finally,  we  find 


S(k,r)  -  8(T)Sn(k,T)'  -  0(-T)SC(k,T) 


(5.31 ) 


'  %  Vk*"^ 
k  " 


*  -itt,  T-Dk  T 


.*  iw^T-Ok  r 
k*ke  K 


•i  I  CkCk«1wk 
k 


,  t  >  t*. 


,  t  <  t'  . 


(5.32) 


5.2  density  of  States 

The- phase  space  volume  Q  corresponding  to  a  variable  $  is  given  by 

a  »  /d*  dir  ■  /dr 


(5.33) 


From  Equations  (4.6)  to  (4.8),  (4.32),  (4.39),  (4.46),  (4.57),  (4.58),  (4.64) 


and  (4.66)  we  see  that 


/dr  -  E 
k 


(5.34) 


By  way  of  illustration,  let  us  assume  an  isotropic  system.  Then  we  have 


I  -  4ir/k2dk  . 
k 


Me  wish  to  convert  the  integral  in  k  space  to  a  frequency  integral: 


4ir/k  dk  -  12)  (w)  du  . 

This  implies  that  the  density  of  states,  2)  (w)  is  given  by 


(5.35) 


(5.36) 


2)  (w)  *  4n  |  k2(w)|  /  I  du/dk  |  .  (5.37) 

Thus  we  can  find  #(«)*  once  we  have  the  appropriate  dispersion  relation.  The 

required  relations  are  given  in  Equations  (4.11),  (4.12),  (4.63),  and  (4.64). 

An  immediate  complication  arises  from  the  fact  that  u  .  and  are  either 

rjk  Ck 

complex  or  pure  imaginary  for  k  $  0,  and  «ak  is  pure  imaginary.  This  means  that 
the  population  of  the  corresponding  states  changes  in  time.  Each  value  of  k 


.  /.  •  .-.v.v.v: 


VV.WL 


corresponds  to  a  unique  value  of  anc*  uak;  so  dispersion  relations 

define  curves  in  the  complex  w  plane.  Let  the  real  part  of  a  frequency  w  be 
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denoted  by  «’  and  the  imaginary  part  by  w".  Then  we  have 


du  _  [  du1  2  ,  du"2  11/2 
dk  *  I  dk  dk 


For  the  A  fields  we  have  the  simple  relation 


<i>"  a  W^a 


(5.38 


(5.39 


Consequently 


-1.5  0.5 

|  kZ  |  dk  »  ^  «a  dw#  . 


If  we  assume  the  existence  of  a  cutoff  frequency  u^,  we  can  normalize 
IMw)  to  obtain 


(5. 40 


2>a(«) 


2  1.5 

um 


(5.41 


For  the  scalar  fields  we  have 


«’k-  /czkz  -  HDk4 


u"  -  ±i0k 
k 


(5.45 


(5.45 


In  Equation  (5.43)  the  +  sign  refers  to  the  C  modes  and  the  -  sign  to  the  17 
modes.  From  Equations  (5.42)  and  (5.43)  we  find 

I-4.I  •  |! 


Using  Equations  (5.38),  (5.42),  and  (5.44)  we  find,  with 


w  ■  Dw/2c 

and  <u  representing  either  u^or 


(5.45) 


f>(w) 


J(W2/  1  -  wz 

/  2  3  4 

/  1  -  3w  ♦  r  w 


(5.46) 


In  the  limit 


we  find 


w  <<  1 


0(w)  - 


(5.47) 


(5.48) 


We  note  from  Equation  (5.42)  that  D  (w' )  has  a  maximum  value  u'  given  by 

m 


#  ^ -  /2c2/d  . 


(5.49) 


This  is  because  the  contribution  to  from  diffusion  causes  the  u’(k)  curve  to 
bend  over  at  large  k,  and  in  fact  to  reach  zero  at  a  critical  value  of  k  given  by 


k"  *  2c/0 
c 


If  0  is  very  small,  we  can  make  the  rough  approximation 

u  »  u'  , 
or 

-  0 


(5.50) 


(5.51) 


(5.52] 


[48] 


a  i  >  A-  k. 


6.  DISCUSSION 


The  approach  taken  in  Section  2  to  handle  the  difficulties  of  dissipation 
is  related  to  the  work  of  Morse,  Feshback,  and  others]^' and  in  fact,  is 
inspired  by  the  earlier  work  of  these  authors.  However,  the  treatment  given 
here  varies  from  earlier  work  not  only  in  that  it  treats  a  more  complicated 
system,  but  in  the  diagonelization  of  H.  Earlier  authors  separated  H  into  a 
conservative  part  that  they  diagonalized  and  a  dissipative  part  that  required 
special  rules  for  incorporation  into  the  theory.  In  Section  4,  the  entire 
Hamiltonian  is  completely  diagonalized,  and  dissipation  is  incorporated  directly 
into  a  formalism  that  is  similar  to  that  of  more  familiar  systems.  This 
approach  is  especially  convenient  for  the  derivation  of  Green’s  functions.  A 
review  of  the  papers  cited  above  and  other  attempts  to  treat  dissipation  in  a 
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Lagrangian  theory  is  given  by  Oekker. 

The  diagonalized  Hamiltonian  for  the  scalar  fields  involves  two  species  of 

functions  in  the  diagonal  or  N  representation.  In  field  theory,  this 

corresponds  to  the  presence  of  a  hidden  variable  that  is  not  explicitly  present 

1 8-21 

in  the  equations  of  motion.  Usually  this  is  a  charge.  Identifying  this 

variable  for  the  scalar  fields  thus  becomes  a  central  problem. 


[49] 
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1 .  INTRODUCTION 


The  primary  use  of  a  field  theory  treatment  for  fluid  flow  will  stem  from 

the  nature  of  the  perturbation  generated.  The  character  of  such  perturbation 

series  is  determined  by  the  relationship  between  the  interaction  terms  and  the 

basis  Hamiltonian,  H  .  In  Part  I,  the  bulk  of  the  formalism  for  a  many-body 

o 

treatment  of  subsonic  Navier-Stokes  fluids  was  developed,  although  the 
discussion  was  confined  to  Hq.  In  the  following  text,  the  nonlinear  terms  in 
the  full  Hamiltonian,  H,  are  considered.  The  discussion  will  be  directed 
towards  flows  with  low  compressibility,  so  that  the  nonlinear  terms  arising  from 
pressure  gradients  will  be  ignored. 

Me  have  succeeded  in  describing  coupling  of  the  scalar  and  vector  potential 
modes.  This  is  done  by  expanding  Hj  in  terms  of  the  diagonalized  fields,  and 
then  constructing  a  linear  transformation  of  those  fields  that  yield  a  new 
rediagonalization,  given  certain  approximations.  The  new  dispersion  relations 
have  a  non-zero  threshold  level  if  the  vector  and  scalar  fields  are  both 
excited.  Also,  the  analysis  of  mode  coupling  due  to  the  nonlinear  terms 
in  the  Hamiltonian  has  been  extended  to  examine  several  similarities  with  the 
literature  on  critical  behavior  and  suggest  a  connection  with  "deterministic 
chaos . " 

2.  MODE  COUPLING 

The  perturbation  Hamiltonians,  H_  and  H...  ,  are  now  included  along  with  H 

I  Nl  O 

to  form  the  total  Hamiltonian  H: 

H  »  Hq  ♦  Hj  +  Hnl  .  (2.1) 

The  objective  is  to  diagonalize  H.  The  Hamiltonian  density,  ,s  easily 
obtained  from  £j,  given  in  Equations  (2.8)  to  (2.12)  of  Part  I.  Me  have 
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M j  *  Jf]  * 


(2.2) 


Jfy  •  2  ly0+  *  *yo)(V#  -  V*)« (VxA  -  VxA) 

Sc 


(2.3) 


»  ~E-^(y0*  *  ^r0)(VxA  -  VxA )  •  ( VxA  -  VxA)  . 
Sc 


(2.4) 


Equations  (2.3)  and  (2.4)  are  expressed  in  terns  of  the  time  derivative  of  the 
generalized  coordinates  rather  than  conjugate  momenta  for  convenience.  An 
alternative  form  for  and  valid  in  the  limit  of  constant  vorticity  is 
discussed  in  Appendix  E. 

It  will  prove  very  convenient  to  use  a  compact  notation  in  the  remaining 
discussion.  The  Fourier-transformed  variables  defined  in  Section  4  of  Part  I 
are  now  expressed  in  the  shorthand  form 


na0no eiq*r  - 

q  nq  q 


(2.5) 


*a+  -ij)«r  + 

n  s„e  y.  -*  n_  , 

q  nq  A  q 


(2.6) 


Vak 


*  _ik*P 


(2.7) 


(2.8) 


The  use  of  the  symbol  T  is  not  meant  to  convey  a  change  from  field  variables  to 
operators,  as  in  quantum  mechanics;  rather  it  is  intended  to  emphasize  that  the 
compressed  notation  includes  y^.  Using  the  transformed  version  of  the  field 


w.v.v.v v.v.w.v  v. -•.*>.>  v.-  .-Va-.-v.-v  v  v.W 


Krnr.'rrv  \n  i  yw  tTi7.im-i.I-t..-. 


variables  given  in  Part  I,  Equations  (4.2)  through  (4.7)  and  Equations  (4.47) 
through  (4.50),  the  following  expressions  are  obtained: 


VxA  -  VxA  ■  i£c/2/p«akkx(a£  ♦  a^)  , 


(2.9) 


Y  7=^  IWn«(nI  *  O  *  “ra^I  " 

°  0  °q  /ptT  ^  d  d  vd  q  q 


(2.10) 


w  -  »*  *  ir,,!  ~  q«J  ♦  c,  -  -lj  -  V  . 


q 


(2.11) 


Y2u 


*  r  + 

7xAro  ♦  rQ7xA  -  crj  —  kx(aR  ♦  a|t)  , 


(2.12) 


Throughout  the  discussion,  frequent  use  will  be  wade  of  the  identity 

/d3r  el(p±q±k)*p  *  «(p  ±  q  ±  k)  .  (2.13) 


We  will  also  take  liberties  with  our  summations,  noting  that  they  go  from  -N  to 
N  with  N  large.  It  will  be  assumed  that 


I  - 1  - 1  . 

K  C  k.q 


N  »  q 


(2.14) 


i.e.,  we  assume  that  the  dominant  interactions  occur  at  long  wavelengths.  We 
begin  with  H^. 

H,  .  JVr  Jt, 

■  *r"  5  lo,,*k  ‘  *f)"’[“Jq*k(Vk  -  Vk>  *  “fq.k'^.k  -  Vk)] 
«;  *  <q  '  "J  -  V/U*k  • 


(2.15) 


Me  also  have 


H2  ■/' 


Jrjt. 


.-fs 


-  1.)  ♦  -  CJJ 


2  u-  cm  j -  l"W"-  V  Cq  -  ‘q' 

16cz  kq  ‘ak  Ypa~  w  q  *  q 


Itkxa^* jx«j6(k  ♦  j  +  q)  +  kxa£»jxaj5(j  -  k  +  q) 


+  kxe^jxepfk  -  j  +  q)  +  kxa£*jxa^6(j  +  k  -  q)] 


H2  reduces  to 


h2  »  -=  U  — <n  -  n!)  ♦  «c  (fq  -  Cl)l 

8/p  kq  Ugj^/wq  nq  q  q  q  q 


j[k*(k  ♦  qje^e. _ -  2k*  (k  ♦  q)*^**  ♦  k*(k  ♦  q)e*a++ 


fk*e _ (k  ♦  q)*a.  -  2k*aJ  (k  ♦  q).eR  +  k*a£  (k  +  q) 

1  k+q  M  ^ 


Finally,  consider  HM,  f  which  is  given  by 


hnl  ■  /' 


d3p  .% 


^)^S=la^  *  “fqlc° 


Wl,* 


♦  *  UD(nO  " 

9  9  9  9  P  P  P 


♦1  »  -  cj)  -  v_(n_  -  C_) 

9  9  9  9  P  P  P 


(2.28) 


(2.29) 


(2.30) 


and  *  “JWkq  '  {P«t)  ‘  w— C—  ’  f— >  •  (2-31) 

K  M  KM  KM  p+q  p+q  p4,q 

In  Equations  (2.19)  to  (2.21)  and  (2.29)  to  (2.31),  as  Mali  as  in  the 
following  material,  the  transformed  Hamiltonians  will  only  be  equivalent  to  the 
original  Hamiltonians  if  the  respective  u's  and  v's  are  equal  to  one. 
Unfortunately,  if  the  transformations  are  to  be  canonical,  the  transformation 


functions  must  be  chosen  to  satisfy 


u’-vi-l 

*  k 


u2  -  v2  -  1 

q  q 


(2.32) 


(2.33) 


(2.34) 


These  requirements  overdetermine  the  u's  and  v's.  A  compromise  may  be  found 


in  which  (letting  j  represent  k,  p,  or  q) 


where 


“j  ■  "j  *  *j 


vj  ■  "j  •  ej 


N.  »  1 


(2.35) 


(2.36) 


(2.371 


t 


I 


The  error  for  finite  N  is  equal  to  1/16Nj  if  the  eigenvalues  of  the  diagonalized 
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terms  in  the  respective  Hamiltonians  are  reduced  by  Nj  for  each  j.  In  the 
remainder  of  this  treatment,  difficulties  with  proper  specifications  for  u^  and 
Vj  will  be  ignored. 

•  * 

Using  Equations  (2.29)  and  (2.31)  H,..  ngy  be  rewritten  as 


hnl  ’  ^  *  *Wc<i '  pMp  *  "> 


.  h3  .  h4  . 


(2.40) 


Similarly,  H2  My  be  expressed  as 

H,  *  H  — 1"  t«  (n  ‘  n-}  +  "ra(fa  "  C-}1  bk+A  k’(k  +  q) 
2  fsTp  kq  «aki^  nq  q  q  Cq  q  q  k+q  k 


*  H5  +  H6  ' 


where  the  shorthand  notation 


k-(k  <•  ■  »’(k  *  1)/%^ 


is  used. 


is  simplified  through  the  change 


-ive 

- 2-  n 

e/T  * 


s  kx(%  - [Vk0; 


q+kcq  ~  ivq+kdq+kdq  ^“ak 


(2.41) 


(2.42) 


(2.43) 


Depending  on  the  eean  amplitude  of  the  A  and  +  fields  it  Mill  be  convenient 


to  group  the  various  terns  in  H  as  either 

H  -  (Hn  ♦  H3  ♦  Hg)  ♦  (Hf  ♦  Ha  *  Hg)  ♦  (Ha  +  H1 ) 


or  as 


(2. 45 


H  H,)  .  (f2Hf  .  Hml)  *  (Ha  ♦  H2>  , 


where  f1  and  ,f2  are  fractions  such  that 


f1  +  f2  *  1 


(2.46 


(2.47 


and  we  define 


H  ■  T  (i) 

n  q  nqqq 

Hc  ■  l  V&  • 


(2.48 


(2.49 


Consider  the  grouping  given  by  Equation  (2.45).  We  write  this  as 

H  «  H*  +  Hr  +  Hs  • 


(2.50 


We  have 


H0  -  k*k 


icy 

kx(v  -  •i'^v4kc,  -  <2-5' 

sic 


£[“ak*k\  -  1,To'*k  '  •l)-Il(Vk"cqlc  ‘  ,Vk"dqk>/B,l 


(2.52 


cy 


Hr  *  *t*Wq',q  *  17=  "m1",  -  "r>  l  k<k  *  q)6£.qV" 


n^J-  IN'  'q  '5'  i  '  '•,"k*q'V”«k 


icy 

— ~ —  u (n_  -  n+)T  p(p  +  q)4+  _♦_/ v'ZTu  ] 

^e/Z^  *»  q  q  p  M  P+«l  P  P  P+q 


^Vj"q  *  V»(nq  *  k(k  *  q,W“«k 

H  ”  ^ 


"  2Vo(nq  "  £  P(P  +  q)  V/VW  ' 


and 


H*  ’  *  ^=*Cq(<q  ‘  k(k  *  W*qV“.k 


icy. 


16/pu^ 


“Cq(Cq  ‘  P<P  *  "'W^] 


*["<qCJCq  *  aJro(Cq  '  ^  k(k  *  ’'W“* 


Kl p<p * ,)nw ]  • 


(2.53) 


(2.54) 


(2.55) 


(2.56) 


2.2  Emery  Transformations 

Decoupled  fields  are  now  obtained  by  applying  canonical  transformations  to 
Equations  (2.48),  (2.50),  and  (2.52)  of  a  type  originally  due  to  Emery. 


is  transformed  by  means  of 


■  Tq  («  .n  .  -  iv  .  n.  .)/w  , 

k  q  q+k  cqk  q+k  dqk  q 


(2.57) 


[61] 


■  klL  « k_ 


.  .»». 


V  ■ il,k2  *  -<k.  , 


(2.58) 


%  "  *k  *  "kV^ak  ' 


(2.59) 


■J  ■  SI  *  ’VoW  • 


(2.60) 


•j  •  *  1fkroY“ak  • 


(2.61) 


In  Equation  (2.54)  we  make  use  of  Equation  (2.52)  and  the  fact  that 


j£  *  -J  , 
*  k 


(2.62) 


as  can  be  seen  from  Equations  (2.47)  and  (2.48).  Using  Equations  (2.55)  through 


(2.57)  in  the  relation 


Hd  ■  £hk*k*k  *  'Vo' v  *j)jk  ] 


(2.63) 


we  find,  remembering  that 


I  V-  "  ^ 

k  k  k  k 


(2.64) 


(Since  the  sums  go  from  -N  to  N),  and  that  =  -1 : 


H*  ■  j[“.k< 


\  -  fkJk/uak 


(2.65) 


This  shows  that  the  normal  modes  corresponding  to  the  vector  potential  in  the 


excited  medium  are  dressed. 


Hr  is  transformed  in  a  similar  manner.  We  rewrite  Equation  (2.53)  as 


Hr  -  I[“nqnjnq  ♦  gqT0(hq  -  nl)[k(k  ♦  qjJ^  -  ip(p  ♦  q)^]] 


(2.66) 


and  make  the  transformation 


n_  ■  r  ♦  9-7  fk(k  ♦  q)J.  -  ip(p  ♦  q)J  ]/u 
q  q  q  o  bq  <Pq  r»q 


(2.67) 


h 


F 


(2.68) 


„t  s  r+  -  g*yo[k(k  +  -  ip(p  +  q)J^]  /w*q 


n+  ■  r+  +  g*YQ[k(k  ♦  qJJ^.  -  ip(p  +  q)J.  J  /«* 
q  q  q  m 


From  Equations  (2.49),  (2. SO),  and  (2.61),  it  follows  that 


Jbq“  \ 


J+  J  .  . 

♦q  - 


Furthermore,  since  g  -  u  , 

q  nq' 


g*/«*  »  g  /u 

V  nq  V  nq 


Using  Equations  (2.61)  through  (2.67)  we  find 


hp  -  Ilvj'-q  +  a5Ck<k  +  q)V  "  1P(P  +  q,V3 


[k(k  ♦  q)^  -  ip(P  ♦  j.  . 


Clearly,  we  may  make  a  like  transformation  for  Hs#  To  wit: 


Hs  *  ^[wCqfJfq  +  flqyoKq  *  “<P<P  + 


Cq  «  sq  ♦  g*r0[k(k  ♦  q)^  -  ip(p  ♦  qN^/% 


Cj  -  *J  *  VoCk(k  +  q,Jbq  '  ip(p  *  q,V  /WCq  ' 


cj  »  s+  +  gqro[k(k  +  qN+q  -  ip(p  ♦  q)J^]  /u>*q  ; 


wv: 


rJWi 


(2.69) 


(2.70) 


(2.71) 


(2.72) 


(2.73) 


(2.74) 


(2.75) 


(2.76) 


leading  to  the  result 


H 

s 


(|>  s*s 

Cd  q  q 


-  g  [k(k  +  q)J 

q  »q 


-  ip(p  ♦  q)J  ] 

<pq 


[k(k  ♦  qJJ^  -  ip(p  ♦  •  (2.77) 

3.  COLLECTIVE  EXCITATIONS 

The  form  of  the  Hamiltonians  H  ,  H  ,  and  H  is  very  similar  to  that  of 

fir  s 

3-9 

electron-phonon  coupling.  This  suggests  that  canonical  transformations  can 

be  made  to  new  Hamiltonians  H',  H',  and  H* ,  which  contain  terms  that  are  only 

p  r  s 

quadratic  and  fourth  order  in  the  fields,  with  no  terms  containing  both  vector 

and  scalar  fields,  and  that  the  new  Hamiltonians  can  be  analyzed  along  the  lines 

8—18 

of  the  BCS  theory  for  superconductivity.  We  now  pursue  each  of  the  desired 
transformations  in  turn.  H  is  given  by  Equation  (2.49): 


H»  -  l  Kk«kV  -  ,Vkkx'*k  '  '  <VkdJ.kV/“q]  •  (3',> 


where 


v  c/*^Sk 


(3.2) 


Equation  (3.1)  may  be  written  as 


»  H  + 
fi  fiO  fil 


(3.3) 


The  transformed  Hamiltonian,  H^,  may  be  written  as 


Hi  *  V  H«  • 
fi  fio  Pi 


(3.4) 


with 


H  »  lu  .  a*a 
fio  £  ak  k  k 


(3.5) 


[64] 


and,  (see  Appendix  F)  again  recalling  that  y*  *  -1, 


Hii 


v 

•k  k 


I 

k 


yr  Mq.k  wg 1  -keo^ccgca 1  -kca '  ,  uq.k  AVl'-k11,' 

«'  "“c,  *  “c,-k»2  *  “V’“5  [(“d,  *  “dq-k»2  *  “Vl"q 


(3.6 


The  meaning  of  the  frequencies  u  .  u  .  . ,  u  .  ,  and  <•>.  .  ,  is  discussed  in 

cq  cq-k'  dq  dq-k 

Appendix  F,  as  are  the  frequencies  u^,  «dq#  *nd  w$q-k*  t0  be  *ouncl 


below. 


Hp  is  given  by  Equations  (2.50)  and  (2.51): 


Hr-I 


jvfo  *  v,<",  -  & 


[l  k(k  .  <)B^V".k  '  2P(P  *  "'♦U*P/VSVq  ] !  '  (3-7 


with 


g„  *  cu  /s/pcT  . 
9  nq  q 


(3.8! 


Me  write  Equation  (3.7)  as 


H  •  H  _  ♦  H  ,  , 

r  rO  rl  ' 


(3.9! 


and  obtain  a  transformed  Hamiltonian  given  by  (see  Appendix  F) 


h; 


H  _  ♦  H’,  , 

rO  rl 


(3.10] 


where 


rO 


r*  f 

l  <*>  n_n„ 

q  W  q  q 


185] 


(3.11 


and 


2  l  flqVjg,P(P  *  q,P’,P'  *  ’'Wp*P'VP'/[(%  '  %-q>  -  "JVw 

-  a  k(k  ♦  q)k'(k'  ♦  ql^Vk'^k^'^bk  -  W2  *  “«•“«} ■  (3-12) 

Me  have  used  a  -  sign  before  the  second  tern  in  Equation  (3.12)  and  in  front  of  w‘ 

within  that  term  as  a  reminder  that  the  frequencies  u. .  and  <*.  are  purely 

OK  bic-q 

imaginary. 


The  expressions  for  H  and  H'  are  similar  to  those  for  H„  and  HI: 

*  s  r  r 

h  «  h  +  h  •  T  u  (+(  +  y  g  ((  -  C+) 
s  so  si  q  l  Cq  q  q  o  q  q  q 


[l  k(k  ♦  '  I*'’  *  ^,VVtVW] |  .  (3-'31 


where 


9q  *  CUCq/8y%  ; 


(3.14) 


h;  *  H,o  ♦  h;i  ' 


(3.15) 


where 


JVic«  ! 


(3.16) 


and 


Hi>  •  - 


( 


g. P(P  *  P.PMP'  *  pt^vJ'-pV'V  ■  W*  *  “CP]V, 


P+q 


*  g-k(k  *  <"k'(k'  *  -  v-,)2  * 


We  now  make  the  definitions: 


mk  ♦  p)b^apk 


bl°a  %  ,  ,  „ 

M  s  k  r//.«  _  /.«  ^  ..2  i..2 


t(v-v-,)  *  %Kk 


0,%  r _ — *  q>BUv 

^  ^  1  k  rfu  _  \2  . 


9  9  *  C(ubk  -  W‘  +  «fqi“ak 


4ck  ■  r**2i  i 


u 


C*  c 


q+k  a+ka 


0  t(“ca  -  2  *  “2k4]«2 


cq  cq-k 


*dk  •  kki  i 


Vq+k  do+kdQ 


q  t(Wdq  '  Wdq-k>2  +  ^ q 


Aiq 


kk  I 


p(p  ♦  q )<p+  rf> 


2*q  q  -  ?  , 

P  “  w.  )  -  w :  ]u  u 

♦p  4>p-q  CqJ  p  p+q 


(3 
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ai 


1*2 
± g  u 
2yq  q 


P  [(«, 


♦P 


p(p  * 

_  Ptq  p _ 

“  W$p-q^  “  wcqjwp“p+q 


(3.23) 


We  now  -interchange  the  pr-iaed  and  unprimed  variables  in  the  double  summations  of 
Equations  (3.6),  (3.12),  and  (3.17),  sue  over  the  primed  variables,  and  use  the 
definitions  (3.18)  to  (3.23)  to  write 

H  *  Ha  +  Hx  '  (3.24) 


where 


"a  ■  l  {“,kVk  -  “*>  ]  ' 


(3.25) 


and 


h*  *  j  )[%❖«.  *  vicq  *  i  [vkci-kV<*  -  'v^AVdk] 


-  I*+  ♦  (A1  +  A"  ) 

P  P-q  P 


(3.26) 


Equations  (3.25)  and  (3.26)  are  now  used  to  find  the  normal  modes  of  the 
coupled  system,  given  the  averages  assumed  by  the  definitions  (3.18)  to  (3.23). 
To  simplify  the  problem,  some  sort  of  approximation  is  needed.  Let  us  choose  to 
replace  the  inner  summations  in  Equations  (3.25)  and  (3.26)  with  typical  values, 
which  we  will  designate  with  capital  letters.  For  example,  we  let 

I-P 

P 

Then  new  canonical  transformations  can  be  made  to  obtain  the  dominant  modes  of 
the  system.  We  first  note  from  the  definitions  (3.18)  to  (3.23)  that  the 
various  As  have  resonant  denominators  at  small  wavenumbers,  so  their  main 


K 


contribution  should  occur  at  long  wavelengths.  Let  us  begin  with  Ha<  With  the 


transformations 


B.k  ■  Bk  *  Vq 


•k  ■  "k  -  N-q 

the  average  value  of  the  second  term  in  may  be  expressed  as 

*  *Sq  ><B!kB.k  -  BV-k>  • 

In  the  limit  as  Q  -  0  and  the  As  remain  finite,  we  find 

H«  ■  -  2(4bo  *  1S0'(*k  -  •i)(*k  -  y]  • 


(3.27) 


(3.28) 


(3.29) 


Again  using  a  standard  procedure  ,  is  diagonalized  with  the  transformations 


where 


“k  *  uk*k  "  v-*-  ' 

k  k 


“k  ’  “k*k  ‘  Y* 


2  2 
u*  -  v  -  1 

*  k 


V  vk  ~  1  • 


(3.30) 


(3.31) 


The  eigenfrequency  for  the  diagonalized  Hamiltonian,  w^,  is  given  by 


where 


wok  "  wok  +  w1k  * 


2  2  2 
wok*wak+w1k  ' 


(3.32) 


(3.33) 


(3.34) 


I 


R 


E* 


Wjfc  now  is  complex  and  is  no  longer  a  simple  quadratic  in  k. 


The  expression  for  the  scalar  Hamiltonian  is  considerably  more  complicated 
than  for  the  vector  Hamiltonian.  Consequently,  it  is  helpful  to  make  a  number  of 
simplifications.  First  we  interchange  the  p  and  q  indices  in  the  last  term  in 
H  .  -We  then  make  the  following  -r,  >roxi mat ions: 

•  let  p,k  0. 

*  *c  *  *d,  -  ‘o  - 

•  Keep  only  terms  with  the  same  value  of  q  in  the  last  ($)  term. 

2 

*  Drop  Dq  compared  to  cq  for  small  q. 

We  then  make  the  definitions 

%  -  W(Adq  +  AJq  >  '  (3*36) 


and 


and  find 


u  ■  lim 
q  q,R-«o 


2<’(u.>  *  V 


v  w 

q  q 


"lq  *  2<"q  "  V  ' 


u  ■  U  ♦  2u  ♦  2A.  ; 

oq  q  q  + 


HX  *  ^“oq(r|inq  *  *  VqCq  *  Vq11' 


(3.37) 


(3.38) 


(3.39) 


Proceeding  as  for  o>a,  we  find  the  eigenfrequencies  of  H  to  be  given  by 


(3.40) 


The  development  given  above  has  strong  parallels  with  the  Bogoliubov  theory 
of  superfluidity8'9, 15-19  and  the  BCS  theory  of  superconductivity.  However, 
the  similarity  in  form  of  Equations  (3.6),  (3.12),  (3.17),  (3.25),  and  (3.26) 
to  the  Hamiltonians  found  in  those  theories  is  not  sufficient  to  indicate 
similar  critical  behavior.  In  the  BCS  theory,  the  Poisson  brackets  are  replaced 
with  (Fermion)  anticommutators,  leading  to  a  dispersion  relation  of  the  fora 

w2*  u2  +  <d2  ,  (3.42) 

O  1 

where  is  a  constant.  This  produces  a  gap  in  the  dispersion  relation  at  low 
wavevectors.  In  the  Bogoliubov  theory,  uQ  scales  as  k2,  with  again  a 
constant.  This  leads  to  a  different  k  dependence  for  the  eigenfrequency  as  k 
increases.  A  dip  in  the  dispersion  curve  can  result  for  appropriate  interaction 
potentials,  creating  an  equivalent  gap  in  the  dispersion  curve. 

In  the  treatment  given  above,  we  have  no  gap  or  dip  in  the  dispersion  curve 
and  hence  no  similar  critical  behavior.  A  gap  will  occur  at  small  k  or  q  if 
the  system  is  bounded,  but  this  boundary  condition  effect  is  different  from  an 
intrinsic  bulk  effect  due  to  many-body  interactions.  Consequently,  no  clear 
conclusion  can  be  reached  regarding  critical  behavior,  as  the  inclusion  of  higher 
order  terms  or  an  examination  of  the  various  As  may  yet  indicate  non-monotonic 
dispersion  curves. 


$3 

4 

•O 
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HAMY-BOOY  TREATMENT  OR  NAVZER-STOKES  FLUIDS 
PART  III,  FURTHER  TOPICS 


[73] 


1.  RULES  FOR  DIAGRAMS 

Rules  for  drawing  diagrams,  summing  them,  and  calculating  the  quantities 

they  represent  can  be  derived  by  examining  the  form  of  the  Green's  function 

given  in  Section  5  of  Part  I  and  following  the  diagram  procedures  used  in 

1  -8 

quantum  field  theory.  The  vector  fields  obey  a  diffusion  equation,  which  is 
very  similar  in  form  to  the  Shroedinger  equation  for  a  nonrelativistic  electron. 
One  problem  that  arises  is  that  the  isolated  disturbances  of  the  solenoidal 
fields  do  not  propagate,  but  simply  diffuse.  The  Green's  functions  for  the 
scalar  fields  are  of  the  same  form  as  those  for  phonons,  which  is  not 

surprising,  since  the  scalar  oscillations  are  essentially  sound  waves. 

* 

Therefore,  the  diagrams  for  the  scalar  fields  should  obey  rules  similar  to  those 
for  phonon  diagrams. 

The  following  rules  are  very  close  to  those  given  by  Shultz:1 

(1)  Draw  all  connected,  topologically  nonequivalent  diagrams  with  2n 
vertices  and  two  external  points,  where  three  lines  meet  at  each 
vertex.  At  least  one  and  as  many  as  three  of  these  lines  may  be  dotted 
lines  with  the  remainder  (zero  to  two)  drawn  as  solid  lines. 

(2)  Write  down  the  contribution  from  each  diagram. 

(3)  Sum  these  contributions. 

The  contributions  to  a  calculated  quantity  corresponding  to  various  diagrams 
are  as  follows: 

(1)  For  the  contributions  to  the  numerator  of  the  two-particle  Green 

function  draw  two  points  labelling  them  r^  and  r£  (at  which  lines  will 
start)  and  two  points  labelling  them  r1  and  r2  (at  which  lines  will  end) 
anywhere  on  the  paper.  Every  diagram  will  ultimately  have  a  line  with 
an  arrow  on  it  leaving  each  point  labelled  with  an  r'  and  a  line 


line  arriving  at  each  point  labelled  with  an  r.  For  the  vacuum 
diagrams,  there  are  no  points  labelled  with  r's  or  r"s. 

(2)  To  get  a  contribution  in  which  interaction  occur  n 2  times,  draw  n2 
points  labelling  them  y1#  y_,....y 

1  c 

(3)  Add  directed  solid  lines  and  dotted  lines  between  points  so  that  each 
"internal  point"  (i.e.,  a  y  or  y')  has  three  lines  either  going  into 

or  out  of  it,  and  each  "external  point"  (i.e.,  an  r  or  r')  has  one  line 
out  of  it  (if  it  is  an  r * )  or  into  it  (if  it  is  an  r). 

•  • 

(4)  For  every  element  yi  y^,  write  G(j(yi|y^). 

(5)  For  every  element  y..  y^.,  write  S(y^  | y^. ) .  Dashed  (phonon) 

lines  need  not  be  directed,  because  (S)  »  S(-k).  It  is  usually 
convenient  to  direct  phonon  lines  just  to  keep  everything  consistent  in 
a  diagram. 

(6)  For  every  element  g  ,  write  po(y_.)  *  lim  G(r.j  -  r2,  -  t) . 

y.j  0  1  t-+o 

(7)  Integrate  over  all  values  from  -«  to  «  of  the  time  and  three  spatial 
coordinates  of  every  internal  point. 

In  practical  flow  fields,  it  will  be  necessary  to  operate  in  coordinate 
space.  However,  since  it  is  usually  much  easier  to  perform  calculations  in 
momentum  space,  the  rules  for  summing  diagrams  in  momentum  space  are  also  given: 

(1)  Draw  all  distinct  structure,  as  in  r-space.  All  vortex  lines  must  be 
directed.  Direct  phonon  lines  for  convenience. 


(2)  Assign  p's  and  q's  to  all  lines  so  that  the  sum  of  the  energy  and 
momenta  entering  a  vertex  equals  the  sum  of  the  energy  and  momenta 


leaving  a  vertex. 


(3)  For  every  element 


include  a  factor  Qq(P) 


(4)  For  every  element  -3— < —  include  a  factor  S(q), 


An  inspection  of  the  nonlinear  terms  in  the  Hamiltonian,  discussed  in  Part 
II,.  shows  that  there  will  be  three  types  of  vertices,  depending  on  whether  the 
number  of  solid  lines  is  0,  1 ,  or  2  as  shown  in  Figure  2.  If  there  are  two 
solid  lines  (Figure  2a),  then  a  vertex  factor  of 


yc  /« 


(1.1) 


should  be  included.  If  there  is  one  solid  line  (Figure  2b),  the  vertex  factor 


» .  He*.  kq  . 

<*2»v  ^ 


(1.2) 


If  there  are  only  dotted  lines  (Figure  2c),  the  vertex  factor  is 


W  »  — - -  p(p  ♦  q)  • 

16/pw 

p+q 


(1.3) 


Since  the  vertices  each  correspond  to  an  odd  number  of  fields,  diagrams  will 
contain  an  even  number  of  vertices.  Assume  that  the  Linked  Cluster  Theorem  holds 


WA'Cv 


v.v.v.-- 


rrrrrin v  v  nr,' vw  vw  v.  v  v  v  ■  ivi-f  n 


and  consider  only  linked  graphs.  The  simplest  two-vertex  graphs  are  shown  in 
Figures  3  and  4.  Apparently,  the  "bubble  graphs"  shown  in  Figure  3  are  not 
allowed  except  when  the  line  joining  the  bubble  to  the  rest  of  the  diagram 
carries  zero  momentum.  Some  4-vertex  diagrams  are  shown  in  Figures  S  and  6. 
Obviously,  the  diagrams  rapidly  become  complicated  as  the  number  of  vertices 
increases. 

Propagators  are  represented  by  uninterrupted  lines;  those  without 
interactions  are  given  by  single  lines,  while  those  with  interactions  are  given 

by  double  lines.  Thus  6  is  given  by  1,  G  by  11,  S  byi  ,  and  S  by  i j 

1  1 1 

2.  SAMPLE  DIAGRAM  CALCULATIONS 

3  9 

Standard  results  can  be  borrowed  from  the  literature  on  diagram  calculations  ' 
in  quantum  many-body  theory,  albeit  some  of  the  diagrams  for  the  Navier-Stokes 
problem  will  have  a  slightly  different  form  than  the  corresponding  condensed 

3 

matter  diagrams.  First,  some  definitions  may  be  adopted  without  alteration. 

A  self-energy  part  is  defined  as  any  diagram  without  external  legs  that  can  be 
inserted  into  a  given  line.  A  self-energy  part  which  cannot  be  broken  into  two 
unconnected  self-energy  parts  by  removing  one  line  is  defined  as  an 
"irreducible"  or  "proper"  self-energy  part.  Let  us  begin  with  Dyson's  equation, 


taking  the  vortex  propagator  as  an  example.  The  procedure  is  to  show  that  the 
sum  of  all  proper  diagrams  can  be  arranged  in  a  geometric  series: 
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Figure  4.  Oyster  Diagraas. 
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Substituting  Equation  (5.7)  of  Part  I  into-Equation  (2.1)  yields 

G  ■  («  -  «ak  +  E)"1  •  (2.2) 

A  siailar  equation  will  hold  for  S.  In  practice  only  a  few  dominant  irreducible 
diagraas  are  suaaed  to  obtain  £.  Fortunately,  a  good  estiaate  for  6  can  be 
obtained  using  Equation  (2.2)  with  only  a  moderately  accurate  expression  for  E. 
Effective  interactions  due  to  sums  over  large  numbers  of  diagraas  can  be 
expressed  in  terms  of  new  types  of  diagrams.  Define  any  diagram  without 
external  legs  which  may  be  appended  to  a  vertex  as  a  susceptance  part,  and  any 
susceptance  part  which  cannot  be  reduced  to  two  simpler  disconnected  susceptance 
parts  by  breaking  a  single  line  as  a  proper  or  irreducible  susceptance  part. 
Exaaples  are  given  in  Figure  7.  Let  the  sum  over  all  proper  parts  be  denoted  by 


Then  in  a  wanner  siailar  to  that  used  to  obtain  Equation  (2.1)  we  find,  if 
bubble  graphs  can  be  ignored. 


Figure  7.  Susceptance  (Unlinked)  Oiagraes. 
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3 .  HEAT  EQUATION 


Me  describe  heat  transport  by  the  equation 
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<7  ♦  +  Oy  V2*?  ■  6 "7  (P  ♦  v*VP)  +  pe 


where 

8'  ■  6/pcp 

0  is  the  coefficient  of  thermal  expansion: 


8 


-  1  (M) 
p  13t;P 


(3.1) 


(3.2) 


(3.3) 


cp. is  the  specific  heat  at  constant  pressure,  and  Dy,  the  coefficient  of  thermal 
diffusivity,  is  given  by 

Dy  *  K/pCp  .  (3.4) 


For  gases  0  -  1/%  while  for  liquids  0  becomes  very  small.  Dy  -jS  usually 
described  by  the  symbol  x,  which  we  are  reserving  for  other  purposes. 

The  symbol 


2 

/ 

where  e  describes  the  generation  of  heat  due  to  friction, 
usually  ignore  it. 


(3.5) 

It  is  small,  and  we 


[87] 


3 . 1  Thermal  Laqranqian 

A  Lagrangian  density  which  describes  Equation  (3.1),  £  ^  is  given  by 
A  “  2(fyo*  '  +  2(V*  -  V^)-(fTT  -  T7T)  ♦  £(?  -  T)2-(72*  -  V2*) 


(3.6) 

(3.7) 

The  third  tern  in  Equation  (3.6)  has  been  inserted  to  avoid  adding  an  extra 

tera  to  the  Navier-Stokes  equation.  The  6'  term  in  Equation  (3.6)  will  add  a 

2 

tera  to  the  Navier-Stokes  equation  of  order  6'0TV“,7  .  We  will  assuae  that  this 

can  be  dropped.  The  last  tera  will  not  contribute  to  Equation  (3.1),  but  will 

1 6 

add  a  tera  of  the  fora 

V7*  70)«  (3,8) 

to  the  Navier-Stokes  equation.  As  in  the  case  of  the  velocity  potential,  the 
siaultaneous  occurrence  in  the  equation  of  motion  of  terms  that  are  of  both  even 
and  odd  order  in  the  time  and  spatial  derivatives  requires  the  use  of  complex 
fields  In  the  Langrangian.  We  will  again  have  relations  for  odd  derivatives 
that  will  require  hypercoaplex  fields  coefficients  to  satisfy: 

7  -  rQf  -  ?rQ  (3.9) 

[88] 


-  tyrf.TT  -  0'TT[*  +  $♦(*-  f*)*(y07*  ♦  v*yq)1 


t-  •  ■  fe’ ! 


*jl(y  t  -  Tro)  -  (7  t0  *  Vo)lfl(*  +  ' 


with 


7  •  T  ♦  T 


.H.J 


s 

V 


and 


VJ  *  YkVT  -  7TYk  .  (3.10) 

If  we  write 

V  -  TkV*  -  7^rk  (3.11) 

and 


P  ■  P(yJ>  *  *YQ)  (3.12) 

and  disregard  terms  that  cancel  upon  variation  of  $  and  ♦  or  are  included  with 
the  velocity  equations,  then  the  lagrangian  L,.  leads  to  the  equations  of  motion 

yQf  ♦  vVT  +  DtV2T  -  /J’T(yoP  +  v»VP)  ♦  |pe  »  0  (3.13) 

and 

-V  -  V-VT  >  0t72T  -  0*T(YoP  +  yVP)  ♦  p>e  -  0  .  (3.14) 

Using  Equations  (3.9)  through  (3.12),  the  addition  of  Equations  (3.13)  and  (3.14) 
leads  to  Equation  (3.1),  as  desired. 


The  momenta  conjugate  to  T  and  T  are  given  respectively  by 


ir_  •  if 

T  2 


(3.15) 


and 


*T  -  -  21  • 


(3.16) 
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The  "Hamiltonian"  density  corresponding  to  £j  is  given  by 


4  *  ^(7*  "  *J>).(i7T  "  TVT)  -  ^(T  -  T) 2 ( V2«|»  -  72*)  +  D.J.7TVT 


♦  fl'TTIf  ♦  *  ♦  W  -  7*).(y  7*  *  7*y)]  ♦  h>AT  *  T)e 

o  o  z  p 


(3.17) 


The  Poisson  bracket  relations  are  similar  to  those  for  the  A  fields: 


[T.(r)(  T^r1)]  =  [T.(r),  T^r')]  »  0 


(3.18) 


[T.(r),  Tj (r' )]  »  £d(r  -  r'). 


(3.19) 


We  may  expand  the  fields  T  and  T  in  a  familiar  way: 


-iw_  t 

T  -  I  Tke  u^r) 


(3.20) 


T  -  I  Tke  IK  uR{r)  , 


(3.21  ) 


where  we  again  take 


uk(r)  ■  e 


From  Equations  (3.15)  and  (3.16)  we  also  have 


Uv 

L  k  * 


iu^t-ik.r 


-  \l  v 

k 


■  ifcu.  t*  ik  •  r 
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The  dispersion  relation  obtained  from  Equations  (3.13),  (3.14),  (3.20)  and 
(3.21)  is 


-iu^  +  iv*k  -  0Tk2  »  0'(rQP  *  V*?P)  ♦  pe/T  .  (3.25) 

If  we  drop  the  tere  on  the  right  side  of  Equation  (3.25),  we  obtain  sieply 

■  v*k  +  iOTk2  .  (3.26) 

Me  see  that  temperature  disturbances  are  daeped  by  diffusion  processes  and  are 
carried  by  the  eean  flow;  they  do  not  actually  propagate,  as  pressure  pertur¬ 
bations  do. 

The  relations 

T  -  [T,  H]  (3.27) 

and 

T  -  [T  ,  H]  ,  (3.28) 

together  with  Equations  (3.17),  (3.18),  and  (3.19)  yield 

y  T  +  V*7T  +  D_V2T  -  0'T(y  £  +  vVP)  »  Jvpe  ♦  (T  -  T)W  (3.29) 

O  T  O  4 

and 

Tyo  ♦  V?T  -  0TV2?  +  8'T(yo^  ♦  v-VP)  -  ^vpe  +  (T  -T)W  .  (3.30) 

Subtracting  Equation  (3.30)  froa  Equation  (3.29),  we  again  obtain  Equation  (3.1), 
as  we  should. 

3.2  General  Treatwent  of  Passive  Scalars 

In  Multicomponent  fluids,  we  encounter  equations  of  the  form 


Equations  of  this  type  nay  be  treated  by  Lagrangian  densities  of  the  fora 


Z£  v  -  Yy  Y  -  Yy  V  ♦  (V*  -  Y*)*(YVY  -  Y7Y)  ♦  (Y  -  Y)2(72#  -  V2*) 

TO  O 


-  2DymY  -  [YfQ(Y)  -  YfQ(Y)]  , 

(3.32) 

where 

#  -  Y  +  Y 

(3.33) 

V  -  V  -  jr  , 

(3.34) 

■  yk?Y  -  7Yyk  , 

(3.35) 

and  fQ  is  the  linearized  part  of  f.  The  conjugate  aonenta  corresponding  to  Y 
and  Y. are  respectively 

irY  -  |y  (3.36) 

and 

frY  «  -  ^Y  (3.37) 

The  Haailtonian  density  corresponding  to  Equations  (3.31)  and  (3.32)  is 
2jf  Y  *  (V*  -  7$)*(YVY  -  YVY)  +  (Y  -  Y)2(V2*  -  Y2#)  +  2DYVYVY 

♦  Yf  (Y)  ♦  Yf  (Y)  .  (3.38) 

o  o 


The  Poisson  bracket  relations  for  the  systea  are  of  the  fora  given  in  Equations 
(3.18)  and  (3.19),  with  T  replaced  by  Y.  We  nay  transforn  Y  and  Y  using  the 
relations 

-iu^t+ik'r 


Y  -  I  Y.  e 


(3.39) 


iu^t-ik^r 
Y  -  I  Yke  ^ 
k 


(3.40) 


[92] 


with 


vk  ♦  iDyk 


(3.41) 


The  description  of  the  total  density  p  is  somewhat  different  than  that  for 
the  species  densities  ^  ,  since  there  is  no  diffusion  term.  Nevertheless,  we 
will  still  need  to  construct  a  Lagrangian  density  in  teres  of  adjoint  fields. 
First  we  show  the  results  of  an  attempt  to  formulate  a  description  in  terms  of  a 
single  field.  The  continuity  equation 


may  be  written  as 


p  ♦  ?• (pv)  ■  0 


p  ♦  Vp-(*F  -  V*)  ♦  p(V2*  ♦  V2*) 


(3.42) 


(3.43) 


Equation  (3.43)  can  be  obtained  from  the  Lagrangian  density  £  given  by 


£p  -  ^[PP  ♦  pVp-W  -  **)]  . 


(3.44) 


Me  see  that  for  a  change  the  density  field  is  purely  real.  makes  no  net 
contribution  to  the  velocity  equations  and  the  equation  of  motion  resulting  from 
it  is  unaffected  by  a  variation  of  the  velocity  fields.  Unfortunately,  however, 
the  momentum  conjugate  to  p,  np,  is  given  by 


(3.45) 


Since  we  cannot  formulate  Poisson  bracket  relations  with  field  variables  that 
are  self  conjugate,  we  must  invoke  adjoint  fields,  even  though  there  are  no 
second-order  derivative  terms  in  Equation  (3.42).  We  see  then  that  Equation 
(3.42)  is  a  special  case  of  Equation  (3.31)  with  =  0. 
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4.  APPLICATION  OF  THE  RENORMALIZATION  GROUP 

by  John  Erdei 

In  this  section  we  Mill  present  initiel  results  on  the  applicability  of  the 
renonaalization  group  to  our  Model.  We  note  that  both  the  static  and  dynanic 
fores  of  RGT  are  of  interest.  The  dynaaic  version  can  be  applied  directly  to 
the  results  of  Appendix  E,  while  the  use  of  the  static  renoraalization  group 
requires  the  recasting  of  the  eodel  in  a  fore  sieilar  to  that  of  Landau-Ginzburg 
Theory.  Such  a  technique  would  be  used  to  study  the  steady-state  features  of  a 
given  flow  configuration.  Keeping  in  line  with  the  thrust  of  the  previous  work, 
we  will  use  the  equations  of  Motion  froe  Appendix  E  to  iepleeent  the  dynaeic 
renoraalization  group  procedure. 

Our  intention  is  to  introduce  a  perturbative  expansion  for  the  equations  of 
notion  given  in  Appendix  E.  The  goal  is  to  categorize  the  fluctuation 
integrals  which  contribute  to  the  renoraalization  of  the  viscosity.  Once  the 
fora  of  the  fluctuation  integral  is  identified,  we  can  examine  the  applicability 
of  the  renormalization  group  procedure  and  determine  a  direction  for  future 
work.  Due  to  the  complexity  of  the  model  and  the  number  of  integrals  contri¬ 
buting  to  the  renormalization  of  the  viscosity,  in  this  report  we  will  display 
the  derivation  of  one  of  the  contributing  fluctuation  integrals.  Two  possible 
expositions  of  the  renormalization  group  can  be  considered;  one  involves  the 

averaging  of  fluctuations  over  all  wavelengths.  The  static  form  of  this  proce- 
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dure  is  described  in  AMit,  while  the  dyneaic  version  can  be  found  in  Lovesey. 

The  second  allows  for  the  averaging  of  fluctuation  over  some  defined  length 

scales.  This  systea  is  in  line  with  the  iterative  form  of  RGT  developed  by 
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Wilson,  and  applied  to  the  Navier-Stokes  equation  by  a  variety  of  authors. 


4.1  Equations  of  Motion 

Me  will  use  the  equations  of  motion  (E.20)  to  (E.2S)  to  generate  a  correc¬ 
tion  to  the  physical  parameters  in  the  systems,  and  then  examine  the  results  for 
the  application  of  the  renormalization  group.  The  complexity  of  these  equations 
of  motion  is  indeed  great,  and  for  this  report  we  will  display  a  portion  of  the 
results  of  a  perturbative  approach  to  fluctuations.  As  a  system  undergoes  some 
transition  to  turbulence,  it  becomes  dominated  by  the  formation  of  vortices. 

The  vortices  are  associated  with  the  local  rotation  of  the  fluid.  As  such,  we 
will  take  the  interesting  fields  to  be  the  set  a  in  the  following  exposition. 
Recall  that  these  fields  are  related  to  the  solenoidal  part  of  the  velocity 
field  A.  Me  begin  with  the  set  of  Equations  (E.20)  to  (E.25).  Me  will  intro¬ 
duce  the  notation  that  any  wave  vector  is  now  placed  as  an  argument  for  each 
field,  leaving  subscripts  to  denote  Cartesian  components.  Mith  this  notation, 
we  have: 

a^k)  »  -iw^a^k)  +  f(k)  J  eil|ik1[A|i(q,k)n+(q+k)n(q)  *  Bm(q,k)n(q+k)C+(q) 


♦  (^(q.kJrjfqKtq+k)  +  0B(q,kK+(q)C(q+k)  ] 


(4.1) 


aj(k)  »  iwaka|(k)  +  f(k)  £  «ill>k1  [AB<q,k)n+(q+k)n(q)  ♦  Bm(q,k)ij(q+k)C+(q) 


♦  CB(q,k)j7(q)C{q«*)  +  D|i(q,k)Ct(q)C(q+k)] 


(4.2) 


n(q)  -  -iw^fKq)  ♦  I  f  <*'  ^np^V*' >  '  > >  [An(q.k'  )n(q-k* ) 


♦  Bn(qfk’)C+(q+k)]  . 


(4.3) 


a.7 


wssssm 


.  ,Vv^V/V.-  W 


+  Cn(qfk(K(q+lca)]  , 

C(q)  -  -<«CqC(q>  ♦  If(k,)<P,tkiIat(k,)  *  ^tBp(q»k' )t|*(q+k' ) 

♦  Dr(q,k')C(q)3 
and 

C+(q)  -  1«CqCf(q)  ♦  I  f(k,)ePatkiCat<k,>  *  *J(k' )]CCp(q»k' )i?*(q*k' ) 

+  Dr(q,k")C+(k’-q)]  . 

Me  Hill  reaove  the  time  derivatives  in  these  equations  by  introducing  an 
appropriate  frequency,  so  that,  for  exaaple,  we  can  set  a  >  -iue.  As  a 
technical  convenience,  we  also  introduce  a  randoa  forcing  term  into  each 
equation.  Using  these  steps,  we  write 

*j(k)  »  0* j- (k) F® (k)  ♦  Q*j(k)f(k)  I  €.1|ik1[Am(q,k)n+(q+k)n(q) 

♦  B|R(q,k)n(q+k)Ct(q)  ♦  C^lq.kJqCqKtq+k) 

♦  D^q.kK+fqKtq+k)]  , 

•j(k)  "  0ij<k)Fi+(k)  +  °?j<k>f <k>  I  €i1*k1tA>(q,k)n+(q+k)n(q) 

♦  B^q.kJnCq+kK^q)  ♦  Ca(q,k)i|(qK(q+k) 


♦  D.(q,k)C+(q)f(q+k)] 


r>(q)  *  Qn(q)Fn(q)  +  Gr(q)  I  f(k')e  k'[a  (k1)  -  a  (k')][A  (q,k* )n(q-k' ) 

k«  npq  p  q  q  n 

♦  Bn(q,k')C+(q+k')]  , 


f(q)  -  Ont(q)Fnt(q)  ♦  Gn+(q)  Z  f(k')e  k'[a  <k')  -  a*(k')] 

jj,  npq  p  q  q 

♦  [An(q,k* )n*(q+k* )  ♦  Cn(q,k‘ ){{q+k' )] 


C(q)  *  GC(q)FC(q)  +  GC(q)  X  f(k,>cp#tk;[«t(k,)  -  aj(k*)] 


♦  [Bt(q.k')q+{q+k')  ♦  Ot(q,k‘K(q)] 


+  (q)  -  GCt(q)FC+(q)  ♦  GC+(q)  (Ic'  > €rstki <»t ’  >  “  aJ<k' 

♦  [Cp(q,k' )q(q+k' )  ♦  Dp(q,k' K+tk'-q)]  . 


•  (- ■ 


aij(k|  ■  Tir^s-T  • 


a  ^q*  "  (-iw  ♦  <»  )  ' 

nq 


(4.15 


a’,t(q)  ■  Hi  -  uTT  • 


(4.16 


°{(ql  -  h«t iS“T  • 


svsl»,V.S^iq?»v»yiv>«?«*?vXv 


V 


“r 


and 


GC+(q) 


(iw  -  i«^q) 


(4.18) 


are  the  bare  Green's  functions  for  the  systee.  We  will  be  Interested  In  working 


with  Equation  (4.1).  Here 


u.k  ■  1vr2 


(4.19) 


4.2  Fluctuation  Integrals 


We  will  impose  the  condition 


<F.(k)F.j(k')>  -  a.yjtkWk-k') 


(4.20) 


where  <  >  denotes  an  average,  and,  if  desired,  some  function  of  k  can  be  defined 
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through  Q(k).  We  will  simply  take  Q(k)  to  be  a  constant.  We  will  introduce  an 
iterative  perturbation  solution  to  Equation  (4.7)  by  averaging  over  the 
variables  n  and  {.  To  do  so,  we  introduce  Equations  (4.8)  to  (4.12)  into 
Equation  (4.1). 

It  is  clear  that  even  with  the  simplification  of  only  treating  the  equations 
for  the  variables  a,  the  complexity  of  the  iteration  is  formidable.  For  the 
purposes  of  this  report,  we  will  examine  one  specific  term  in  the  series,  and 
show  its  relationship  to  the  formal  renormalization  of  v.  Under  the  described 
iteration  procedure,  one  can  show  (See  Appendix  G)  that  Equation  (4.7)  can  be 
recast  into  the  form 


<aj(k)>  -  <G*j(k)F*(k)> 

♦  <Q®j0«)€wlB|k1erstks  I  Gn(q+k)Gn+(q)Gr?(q+k) 

tf(k)Am(a,k)3[f(k)Ar(q,k)]G*t(k)]F*(k)>  ♦ 


(4.21) 


I  Qn(q+k)Qr,+  (q)Qn(q+k)f(k)A  (q,k)f(k)A  (q,k) 

q  *  r 

Q*t(k)}F*(k}  ♦  .  (4.22) 

Me  can  therefore  define  a  (partially)  renorealized  Green's  function  in  the  form 

®ij(k)  «  <j<k>  ♦  G?t<k>Wrstklksf(k)Qit<k> 

Xa  (q,k)A  (q,k)G^(q+k)Gr^‘(q)G,?(q+k) .  (4.23) 

q  "  r 

One  of  the  effects  of  the  non-linear  coupling  between  fields  is  to  alter  the 
matrix  structure  of  the  bare  Green's  function,  so  that  it  becomes  non-diagonal. 
The  off-diagonal  elements  indicate  that  the  viscosity  matrix,  which  was 
originally  diagonal  in  the  linear  theory,  becomes  a  non-diagonal  viscosity 
matrix.  This  can  be  seen  since  one  can  invert  Equation  (4.23)  and  identify  the 
viscosity  matrix.  In  the  linear  theory. 


CQ® j <k) ] ”1  -  ♦  i«ak) 


-  -iw^j  -  C»«i j 3k 


(4.24) 

(4.25) 


which  defines  the  viscosity  matrix  v«„.  When  the  (partially)  renormalized 
Green's  function  given  in  Equation  (4.23)  Is  inverted  and  the  renormalized 
viscosity  identified,  the  viscosity  matrix  will  no  longer  be  diagonal,  i.e., 


.a  _i 
[8^(10] 


-ialn  ■  vijk  ■ 


(*.26) 


One  night  associate  the  appearance  of  off-diagonal  elements  in  the  renormalized 
viscosity  with  the  generation  of  a  turbulent  eddy  viscosity  matrix.  The  exact 
nature  of  these  off-diagonal  elements  requires  further  study. 


We  see  froe  Equation  (4.23)  that  a  complete  description  of  the  viscosity 


matrix  requires  the  evaluation  of  the  sum 

I.r(k)  «  lA-(q,k)Ar(q#k)Qn(q+k)Qn+(q)Qn(q+k).  (4.27) 

Mr  q  m  r 

This  sum  typifies  one  of  the  contributions  of  the  i)  and  C  fields  to  the  study  of 
the  properties  of  the  a  fields.  It  mould  be  convenient  for  future  evaluations  to 
replace  this  sum  with  an  integral, 

I  -  f  d3kg(k),  (4.28) 

k  J 

where  g(k)  is  an  appropriately  defined  density  of  states.  Then  the  sun  given  by 
Equation  (4.27)  can  be  written  in  the  form 

Xmr(q)  *  l  d3«<q> 

[(2qw  +  kw»aa(k)  +  (q,  ♦  ♦  q,<*n(q)] 

[«  -  «n(q-Mc)]2 


[ (2qp  ♦  kr)ua(k)  +  (qr  +  k^u^q+k)  ♦  qpun(q)] 


[«  -  «n(q)] 


(4.29) 


or. 


I_r(<l)  *  2f  dx  f  d3qg(q)[(2q  ♦  k  )«  (k)  ♦  (q  +  k  )u  (q+k)+q  u  (q)] 
■r  J  o  J  mna  m  m  rj  mi? 


E (2qr  ♦  kp)u>a(k)  ♦  (qp  ♦  kp)un(q+k)  ♦  qpun(q)] 
t«  -  (1  -  x)«  (q)  -  xu  (q+k)]3 


(4.30) 
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where  we  have  used  the  Feynman  Method  of  folding  denominators 


*1  a“2  a*" 

“l  *2  “•  *n 


Ha,-  +  a2  ♦  ...  an) 
Ha,)  r(a2)  ...  r«xn) 


|  dx, 


dx2-"‘1Vl 


a,-1  a_-1  a  a  -1 

*i  x2  Vi  (1  -  x  -  x  -  ...  -  »n.,) 


[«1*1  '♦  x2*2  *  •"  ♦  Vi  Vi  *  <’  -  X1  -  x2  -  •  Vi".1 


(4.31) 


0  <X  <  I  ;  X  +  X  +...+X  <1. 

i  1  2  n-1 


(4.32) 


At  this  point,  the  requirement  of  a  renormalization  group  calculation  can  be 


based  on  two  points.  If  the  Integral  I>r(q)  suffers  from  infrared  (q  -  0) 
divergences,  then  the  renormalization  group  method  can  be  used  to  discuss  the 
removal  of  divergences.  One  may  also  incorporate  a  procedure  which  allows  for 
the  treatment  of  the  fluctuations  in  the  q  and  C  fields  only  over  a  prescribed 
set  of  length  scales.  Assuming  a  constant  density  of  wavevectors,  the  integral 
I  (q)  indicates  a  q-dependence  of  the  form 

Ml 


q  ,  (4.33) 

which  suggests  that  the  integrals  will  suffer  a  power-law  divergence  for  spatial 
dimensions  below  six.  A  renormalization  group  treatment  which  employs  an 
e-expansion  technique  would  be  valid  only  around  six  dimensions,  and  would  not 
be  reliable  for  a  three-dimensional  system.  Similarly,  one  would  have  to 
determine  if  the  iteration  is  defined  such  that  higher-order-fluctuation 
integrals  do  not  introduce  higher-order  divergences.  If  so,  the  renormalization 
group  procedure  will  f ai 1 . 


The  procedure  nay  be  used,  however,  if  one  assumes  an  upper  and  lower 


cut-off  in  wave  vector.  For  example,  shell  integrations  may  be  used  to  integrate 
over  length  scales  from  a  long  wavelength  on  the  order  of  the  size  of  the 
appropriate  region,  down  to  the  Kolomogorov  scale.  Me  could  then  carry  out 
integrals  of  the  type 


f  A  *  f^*1  ■  <*-34> 

where  the  A's  define  the  wave-vector  shell.  The  shell  integration  could  then  be 
iterated  over  the  shells  (i.e.,  over  values  of  i)  until  the  entire  wave-vector 
region  is  covered,  or  the  effects  of  fluctuations  on  a  specific  scale  can  be 
treated  by  examining  only  one  of  the  shell  integrations.  This  method  would  not 
result  in  divergent  integrals  due  to  wave-vector  dependences.  Since  the 
interactions  between  the  a  fields  and  fluctuations  of  specific  lengths  are  of 
interest,  this  is  the  technique  we  would  pursue. 
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SUMMARY  AND  SUGGESTIONS  FOR  FUTURE  WORK 


Me  now  review  the  developaent  of  the  work  described  in  Parts  I  through  III 

i 

and  explore  the  implications  for  further  study. 

Su— arv 

The  lagrangian  for  corresponding  to  the  Navier-Stokes  equation  was  presented 
in  Part  I,  Section  2,  and  the  basis  for  deriving  it  is  given  in  Appendix  A. 
Equations  of  notion  for  the  linearized  system  were  derived  using  Euler-Lagrange 
relations.  The  Hamiltonian  and  Poisson  bracket  relations  were  derived  from  the 
Lagrangian,  and  equations  of  motion  were  again  obtained.  The  equations  of 
motion  obtained  from  the  Hamiltonian  and  the  Poisson  bracket  relations  were 
identical  to  those  obtained  from  the  Lagrangian,  showing  consistency  in  the 
theory.  The  relationship  between  the  above  equations  and  the  Navier-Stokes 
equation  was  also  demonstrated.  The  entire  developaent  given  in  Section  2  of 
Part  I  depends  upon  the  assumption  that  the  velocity  potentials  may  be  expressed 
in  terms  of  the  sums  of  hypercoaplex  fields  and  their  adjoint  fields.  The 
justification  for  this  assumption  is  discussed  in  Appendix  C.  It  is  not  clear 
whether  the  postulates  given  in  Equations  (2.5)  and  (2.6)  of  Part  I  contain  a 
physical  basis  or  are  merely  a  mathematical  device.  However,  without  these 
postulates,  the  entire  theory  falls  through.  In  Section  3  of  Part  I  the 
symmetries  of  the  Lagrangian  were  examined  to  derive  conservation  laws  and  to 
determine  those  variables  which  are  conserved  with  respect  to  the  basis 
Hamiltonian,  Hq.  The  derivation  of  the  important  stress-energy  tensor  is  given 
in  Appendix  0.  The  form  given  in  Appendix  D  pertains  only  to  the  scalar  fields. 
A  four-vector  formalism  was  derived  to  include  the  solenoidal  fields  in  the 
September  1985  monthly  report,  but  was  not  included  in  this  report  as  it  was 


judged  to  be  unwieldy.  The  symmetry  properties  of  the  Lagrangian  and  conserved 
variables  can  be  used  to  analyze  the  transition  from  laminar  to  turbulent  flow. 

The  Hamiltonian  was  diagonalized  in  Section  4  of  Part  I.  It  was  found  that 
two  different  fields  are  required  to  diagonalize  the  scalar  fields.  This 
suggests  an  analog  with  the  description  of  a  charged  meson  in  quantum  field 
theory  or  perhaps  the  electron-hole  description  of  a  metal.  However,  the 
physical  significance,  if  any,  behind  the  need  for  two  different  diagonalized 
fields  remains  elusive.  One  possibility  is  that  the  fields  describe  sound  waves 
in  which  local  compressions  are  in  phase  with  the  maxima  in  the  local  velocity, 
and  that  the  (  fields  describe  sound  waves  in  which  local  condensations  are  out 
of  phase  with  the  local  velocity. 

Green's  functions  were  derived  in  Section  S  of  Part  I.  The  form  of  the 
propagators  for  the  potential  fields  resemble  phonon  propagators,  while  the 
propagator  for  the  vortex  fields  resembles  that  of  a  nonrelativistic  electron. 
This  suggests  that  there  nay  be  an  analogy  between  fluid  dynamics  and  the 
electron-phonon  problem.  The  development  of  Part  I  contains  all  the  elements 
needed  for  a  field  theory.  Part  I  thus  describes  the  core  work  of  this  report. 
Its  consequences  were  explored  in  Parts  II  and  III. 

Part  II  illustrated  how  the  canonical  transform  procedure  common  to  quantum 
field  theory  can  be  used  to  investigate  mode  coupling  and  reveal  analogies  with 
various  critical  systems  which  are  relatively  well  understood.  The  discussion 
of  Part  II  suggests  that  the  many-body  formalism  can  be  useful  not  only  for 
describing  the  interactions  between  distrubances  in  fluid,  but  can  also  be  used 
as  the  basis  for  other  descriptions  of  turbulence,  such  as  the  "deterministic" 
approaches  currently  in  vogue. 


Part  III  is  essentially  a  preview  of  more  advanced  calculations  which  can  be 
based  on  the  developments  of  Parts  I  and  II.  Rules  for  diagram  calculations 
were  given  in  Section  1  of  Part  III,  and  simple  illustrative  calculations  were 
given  in  Section  2.  The  rules  were  made  by  analog  with  the  rules  for  condensed- 
matter  theory,  and  may  need  to  be  modified  following  a  thorough  review. 

Section  3  of  Part  III  showed  how  the  concepts  of  Part  I  may  be  applied  to 
describe  additional  phonemena,  especially  the  behavior  of  passive  scalars.  The 
heat-budget  equation  was  used  as  an  example.  This  treatment  shows  that  the 
concepts  used  in  this  report  can  be  applied  to  the  solution  of  a  wide  range  of 
partial  differential  equations,  including  but  not  limited  to  other  dissipative 
motion.  Section  4  gave  an. example  of  a  renormalization  group  calculation  which 
was  done  by  John  Erdei.  This  extensive  calculation  illustrates  how  the 
renormalization  group  can  be  used  to  include  multiple  scales  in  a  calculation. 

Conclusions  and  Suggested  Research 

The  formalism  developed  in  this  report  is  internally  consistent,  and  shows 
strong  parallels  with  field  theories  for  other  systems,  as  was  intended.  It 
appears  to  satisfy  the  objective  of  deriving  a  Hamiltonian  for  the  Navier-Stokes 
equation  that  can  be  used  to  carry  out  many-body  calculations  of  turbulent 
flows.  This  work  indicates  that  all  of  the  standard  techniques  of  field  theory, 
including  canonical  tranformations  to  rediagonalize  the  field  variables,  the  use 
of  propagators,  and  diagram  calculations  can  be  brought  to  bear  on  turbulent 
systems  using  this  or  an  allied  formalism,  notwithstanding  that  the  development 
is  entirely  classical,  and  includes  dissipation.  Indeed,  the  approach  taken 
here  is  applicable  to  a  wide  range  of  other  partial  differential  equations. 


A  number  of  loose  ends  remain  in  the  presentation  of  the  basic  theory.  In 
particular,  the  Green's  functions  and  rules  for  diagrams  should  be  thoroughly 
reviewed  and  derived  on  a  rigorous  basis.  In  addition,  the  symmetries  of  the 
Lagrangian  should  be  analyzed  further,  particularly  in  regard  to  local  gauge 
invariance. 

To  merit  serious  attention,  any  theory  oust  give  reasonably  accurate 
calculations  of  experimental  results.  Therefore,  the  present  formalism  should 
be  tested  as  soon  as  possible  before  extensive  work  is  based  upon  it.  Since 
almost  all  calculations  will  involve  Green's  functions,  diagrams,  or  probability 
density  functions  and  partition  functions,  these  quantities  should  be 
rigorously  derived  and  checked  for  internal  consistency.  The  nonlinear 
expansion  terms  in  the  pressure  gradient  term  were  not  studied  in  terms  of  the 
diagonalized  field  variables.  The  contribution  of  these  terms  to  mode  coupling 
should  be  examined. 

Once  the  foundations  of  the  theory  have  been  checked,  calculations  should  be 
made  of  those  aspects  of  turbulence  which  are  well  understood.  Turbulence 
spectra  should  be  predicted,  the  propagation  of  sound  through  a  turbulent  fluid 
should  be  analyzed,  and  expressions  for  the  lower-order  velocity  moments  should 
be  derived.  Finally,  attempts  should  be  made  to  predict  the  transport  of  mass, 
momentum,  and  heat  for  a  turbulent  fluid  in  a  simple  geometry. 
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APPENDIX  A 


EULER- LAGRANGE  EQUATIONS  FOR  FIELDS 


Consider  the  functional 
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Using  Equations  (A. 3)  to  (A. 7),  the  condition 

6L/£a  »  0 

leads  to  the  equation 
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Relation's  similar  to  (A. 4)  to  (A. 7)  hold  for  the  vector  fields,  A.  For  example 
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Continuing  in  a  parallel  manner  to  that  used  to  derive  Equation  (A. 6)  to  (A. 9), 
the  condition  (A. 8)  gives 


U  .  i-M  .  ,x-lL  ,  i_  vx-^4-  .  o  . 
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APPENDIX  B 


OPERATOR  FORMALISM 

As  noted  earlier,  the  solution  to  any  partial  differential  equation  (PDE)  is 
a  field.81  vie  add  that  it  is  a  field  under  scalar  addition  and  eultiplication 
over  the  domain  of  the  independent  variables  of  the  PDE.  In  our  case,  this 
domain  is  {r,t}.  In  passing.  Me  also  note  that  the  field  variables  for  each  PDE 
fore  an  Abelian  group  under  scalar  addition92  and  that  the  Poisson  bracket 
relations  for  these  fields  and  their  conjugate  fields  define  a  Lie  algebra.83 
Of  eore  iaaediate  interest  is  the  fact  that  any  function  of  a  field  variable 
defines  a  vector.  In  particular,  the  elements  of  the  fields  themselves  comprise 
a  vector  through  the  relation 

f(*)  -  ♦  ,  (B. 1 ) 

where  we  now  use  the  letters  ijr,  and  x  to  denote  general  field  variables. 

Now  for  every  vector  field  there  exists  a  dual  vector  field.  We  may  think  of 
vectors  defined  by  Equation  (B.1)  as  column  vectors  and  the  dual  vectors  as  row 

B2 

vectors.  Then  the  dual  vectors  are  formed  by  taking  the  transpose  of  the 
orginal  vectors. 

We  may  take  the  outer  or  tensor  product  of  a  vector  and  its  dual  vector  to 
form  a  tensor  that  can  be  represented  by  a  matrix.  Dirac  invented  a  handy 
notation  to  take  advantage  of  the  above  properties  of  field  variables.04  His 
application  was  confined  to  the  solution  of  the  Wave  equation  for  the 
probability  amplitude  in  quantum  mechanics,  but  we  may  use  it  quite  generally. 

We  denote  a  column  vector  by  the  "kef'  j$>,  and  its  dual  row  vector  by  the  "bra" 
<+|.  Then  the  tensor  product  of  a  column  vector  times  a  row  vector  is  denoted 
by  |t><t|t|.  P*<*t  of  the  utility  of  this  matrix  formalism  is  due  to  the 
isomorphism  that  exists  between  the  properties  of  linear  operators  operating  on 


fields  and  the  properties  of  matrix  multiplication.  This  happy  correspondence 
is  further  strengthened  by  the  fact  that  the  isomorphism  extends  to  the  trans¬ 
formation  properties  of  these  quantities.  However,  if  the  matrix  formation  is 
to  be  truly  useful,  we  must  define  an  inner  product,  which  we  denote  by  <$|tp>, 
that  has  the  following  properties:82 
to  linearity 


<c0  +  <*Hx>  *  c<*|x>  +  d<#|x>  , 


(B.2) 


skew-symmetry 


and  positiveness 


<♦!*>*  «  <*|<t»  , 

«fr|<fr»0,  real,  $  $  0 


(B.3 ) 


(B.4) 


In  Equations  (B.2)  to  (B.4)  we  have  denoted  complex  scalars  by  the  letters  c  and 
d.  If  the  fields  4  and  are  real,  then  the  skew-symmetry  property  collapses  to 
simple  symmetry.  The  properties  defined  by  Equations  (B.2)  and  (B.3)  give  the 
inner  product  the  property  of  bilinearity. 


If  we  can  find  an  operation  that  meets  the  requirements  of  Equations  (B.2) 
to  (B.4)  we  will  have  vector  fields  that  form  a  unitary  space,62  and  will  have  a 
complete  formalism  that  can  be  used  to  calculate  quantities  of  interest.  Many 
linear  operations,  such  as  conditional  probability  density  functions,  satisfy 
Equations  (B.2)  and  (B.4).  However,  it  is  more  difficult  to  satisfy  Equation 
(B.3).  The  requirement  of  Equation  (B.3)  is  met  by  the  definition 

«fr(a)|*(a)>  ■  /da  $*(a)#(a)f (a)  ,  (B.5) 

where  a  is  the  common  domain  of  the  fields  $  and  $  and  f(a)  is  the  steady-state 
probability  density  function  for  the  occurrence  of  given  values  of  a.  Equation 
(B.5)  defines  a  correlation  function,  which  is  a  very  useful  quantity  in 


[114] 


probabilistic  descriptions  of  many-body  systems.  The  utility  of  the  definition 

(B.5)  Mas  first  pointed  out  by  Zwanzig  85-812  ancj  Mori81 3-617  and  has  led  to  the 

si a. non 

development  of  projection  operator  and  memory  function  calculations.  ° 

Actually,  the  definition  (B.S),  aside  from  the  Meighting  factor  f(a)  that 

takes  account  of  the  possibility  expressing  the  fields  +  and  9  over  an  arbitrary 

transformed  domain  for  which  all  values  a  say  not  be  equally  probable,  is  the 

B2 

same  as  a  standard  mathematical  definition  of  the  inner  product. 

We  see  that  we  may  construct  a  matrix  formalism  that  is  completely  parallel 
to  the  well  known  one  of  quantum  mechanics.  This  formalism  is  quite  general  and 
not  at  all  peculiar  to  a  quantum  system.  Indeed,  almost  all  of  the  development 
stemmed  from  the  basic  fact  that  the  solution  to  any  POE  is  a  field.  The  only 
aspect  of  our  development  that  presented  any  serious  restriction  was  the  skew- 
symmetry  requirement  for  the  inner  product. 

Since  they  are  vectors,  |+>  and  <$|  satisfy  the  relations 


|$>  ♦  |«|r>  »  |<Jr>  ♦  |$>  ,  (B. 6) 

«H  ♦  <«fr|  ■  ♦  «*|  ,  (B. 7 ) 

i+>  ♦  (|l»>  ♦  lx>)  ■  (i*>  +  |*>)  +  |x>  ,  (B.8) 

<♦1  ♦  («fr|  ♦  <x|)  *  (<♦!  ♦<♦!)♦  <x|  r  (B.9) 

c(|$>  ♦  |*>)  -  c !♦>  ♦  c Idr>  ,  (B.  10) 

c(«*|  +<♦!)*  c«>|  +  c<*|  ,  (B. 1 1 ) 

and 

i|$>  »  | q»  ,  i«fr|  ■  <$| 


Furthermore,  the  dual  of  c|*>  is  c*<$j.  we  write  this  as 


(B. 12) 


Let  us  denote  the  linear  tensor  operators  |$><$|,  |$><$|,  and  |x><x| 
respectively  by  *,  T,  and  x.  Since  the  transformation  properties  of  these 
operators  are  isomorphic  to  those  of  matrices,  we  also  have  the  relations 


c(d#}  »  (cd)* 

(c  ♦  d)*  »  c*  ♦  d* 
c(*  ♦  f)  *  c*  +  cf 
♦(*X)  »  («)X 

(*  +  *)x  *  *x  +  Tx 
*<*  *  x)  *  w  ♦  *X 
(c*)f  a  c(*f) 
and 

♦(cf)  »  c(«) 


(B. 15] 
(B. 16) 
(B.  17 
(B.  18 
(B.  19 
(B.  20 

(B.  21 


Furthermore,  there  will  exist  an  identity  operator,  I,  such  that 

!♦  *  ♦  »  *1  (B.22; 

The  combination  of  Equations  (8.6)  through  (B.13)  and  (B.14)  through  (B.22) 

a  a  do i  poo 

leads  to  the  usual  linear  relationships.  '  '  The  operation  of  a  bilinear 

operator  on  a  vector  will  yield  a  new  vector  of  the  same  type, 

♦  |<F>  *  |x>,  <*|*  3  <x|  •  (B. 23 ] 

This  is  often  taken  as  the  definition  of  a  bilinear  operator . 84 » 821  * 822 
Me  also  have 


♦(!♦>  +  |x>) 

a  ♦(»!»>  ♦ 

♦  lx>  , 

(B. 

24 

(<*l  +  <xj# 

a  + 

<x|*  , 

(B. 

25 

*(c|<F>) 

a  c*  |<fr> 

9 

(B. 

26 

(<«j»|c)* 

a  <#|4c 

9 

(B. 

27 

(♦  ♦  T)|x> 

*  »jx>  ♦ 

t|x>  , 

(B. 

,28 

<x|(*  +  T) 

=  <x|*  + 

<x|f  , 

(B. 

29 

(B.  2 


and 


(♦f)|x>  *  *(T|x»  , 
<x|(**)  *  (<x|*)f  . 


(B.30) 


5WS 


However,  we  do  not  necessarily  have  coMiutation  of  operators; 

W  t  tt 


(B.31) 


(B.32) 


Note  that  all  of  the  properties  (B.6)  through  (B.32)  are  quite  general,  and 
beginning  with  equations  like  Equation  (B.1).  can  be  derived  for  any  field.  The 
critical  properties  (B.2)  to  (B.4)  depend  on  the  definition  chosen  for  the  inner 
product . 


We  now  define  the  length  |*{  Of  a  vector  |*>  or  <*|  by82 

1*1  ■  <*!*>* 


(B.33) 


If  all  vectors  have  a  finite  length,  then  Equations  (B.2)  through  (B.33)  define 
a  Hilbert  space.  If  in  addition,  our  definition  of  the  inner  product  is  such 
that  the  rules  (B.2)  to  (B.4)  apply,  we  eust  have:82 

|c*|  *  |c||*|  ,  (B.34) 

|*|  >  0  ,  *  *  0  ,  (B.35) 

|<*|*>|  <  |*||*|  ,  (B.36) 

I*  ♦  *|  <  1*1  +  |*|  ,  (B.37) 

I*  -  *|  •  0  ,  (B.38) 

|*  -  *|  >  0  ,  *  *  *  ,  (B.39) 


(B.40) 


!*  -  *|  ♦  I*  -  x|  >  I*  -  x| 

Equations  (B.34)  through  (B.40)  define  a  Euclidian  vector  space. 

Froe  Equations  (B. 13)  to  (B.21)  we  see  that  we  must  also  have 

<*|  (c |*> )  ■  c<*|*>  , 

<*|(|*>  +  |x>)  *  <*|*>  ♦  <*|x>  , 

(<*|  ♦  <*)|x>  3  <*(x>  ♦  <*|x>  . 
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(B.41 ) 
(B.42) 
(B.43) 


These  Equations  (B.41)  to  (B.43),  as  well  as  aany  of  the  equations  above  are 
given  as  postulates  by  Fujita,®^!  ^  see  that  they  follow  directly  fro*  the 
properties  of  fields  and  our  definition  (B . 5 ) . 
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APPENDIX  C 


HYPERCOMPLEX  DESCRIPTION  AND  COUPLED  FIELDS 

Me  now  elucidate  the  incongruous  situation  which  occurs  whenever  we  have 
first-order  tine  or  space  derivatives  in  our  equations.  This  occurs,  for 
example,  when  we  juxtapose  the  equations  of  notion  for  the  tft  and  ♦  fields  or 
the  A  and  A  fields,  and  is  most  evident  in  the  disquieting  form  of  Equations 
(2.38)  and  (2.39)  in  Part  I: 

♦  *  yjb  -  ir0  »  (C.i) 

and 

m  »  -  tyyk  .  (C.2) 

We  may  find  some  relief  from  the  unfamiliar  appearance  of  such  equations  by 
noting  that  the  adjoint  fields,  ♦  ,  defined  by 

*  *  **rA  ,  (c.3) 

where  yA  is  a  suitable  hypercomplex  number,  behave  differently  than  the  more 
familiar  fields,  such  as  «|»*.  The  use  of  hypercomplex  numbers  is  common  in 
particle  physics,  where  they  are  used  in  a  four-dimensional  context  to  describe 
degrees  of  freedom  of  the  system,  such  as  spin,  that  are  not  explicitly  manifest 
in  the  equations  of  motions.  The  physics  underlying  the  use  of  hypercomplex 
numbers  in  our  equations  is  not  yet  clear,  but  a  hypercomplex  description  gives 
us  extra  freedom  to  achieve  a  fully  self-consistent  system  of  equations.  Me 
first  give  a  brief  review  of  hypercomplex  numbers,  and  then  illustrate  how  they 
may  be  used  to  deal  with  anomalies  in  the  equations  of  motion. 

Hypercomplex  numbers  are  an  extension  of  complex  numbers  and  are  often 
represented  by  matrices,  as  their  multiplication  table  can  be  illustrated  by 

1120] 


V  v  “>  v  % 


k  I'ft  IV,  I,  ll. 


,4  ||  <‘afc 


Multiplication  of  tha  appropriate  matrices.  In  our  context,  they  are  under¬ 
stood  to  operate  in  a  four-space  apart  from  the  space- time  space  described 
above.  Me  are  especially  interested  in  four  such  numbers, 


V  V  V  y3  "  V  T  \  ' 


which  we  will  now  associate  with  the  time  and  space  derivative 

a,  ■  v  ax  • 


as  follows: 


V  ax  -  W  Vx  • 


(C.4) 


(C.5) 


(C.  6) 


We  have  a  wide  degree  of  latitude  in  choosing  matrices  to  describe  the  numbers 
y ;  the  choice  will  rest  on  those  aspects  of  the  system  we  wish  to  highlight.  4 


The  Pauli-Oirac  representation: 


k  a 


where  the  o..  are  the  Pauli  spin  matrices. 


(C.7) 


(C.  8) 


i  L  1 


i  1 

,  a  *  ,  O  a 

2  i  3  1 


(C.9) 


is  the  most  commonly  used  system  in  particle  physics,  and  emphasizes  parity 
conservation.  The  Weyl  representation: 


(C.10) 


is  also  common  to  particle  physics,  especially  when  considering  the  conservation 
of  chirality. 
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liarflb 


[?: 


We  will  define  a  representation  that  will  emphasize  symmetry  properties 


under  time  reversal: 


We  find 


and 


1 

#  y  ■ 

a 

1 

k 

a 

Vk 


■  Vo 


(C.11) 


(C. 12) 


(C. 13) 


(C.  14) 


We  now  postulate  a  pair  of  fields  $  and  x,  coupled  by  the  first-order 

C2  C3 

hypercoaplex  equations.  ' 

<Toat  *  C1,k3x>*  *  *°<Vt  -  cVx>Vx<»  -  *>  •  «c-,s> 


and 


(To*t  *  clrkax>*  ■  i°<Vt  *  'Vx'W*  ■  x)  •  lc-,6) 


Operating  on  the  first  of  these  equations  by  yQ3t  ♦  cy^  and  on  the  second  by 


*Vt  *  =Vx 


<®t  -  -=*ax>*  ■  i°(at  -  WA  *  c2^»Vx<*  •  *1 


(C.17) 


and 


<8t  •  cZax>*  ■  4(a»  *  !VoVi  *  c‘aXv*  •  *>  •  <c-is» 


2.2, 
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Subtracting  Equation  (C.18)  froa  Equation  (C. 17),  leads  to 


<*t  -  cV)  (♦  -  X)  •  -  owAW*  -  x> 


■  W,"«»  -  *>  • 


(C. 19) 


If  we  now  sake  the  connection 


♦  * 


(4) 


(C.20) 


we  retrieve  our  equation  of  Motion  for  t  corresponding  to  L  : 


(at  "  «2*>  *  PVxVx* 


Froa  Equation  (C.21)  we  find  the  corresponding  Equation  for  #  ; 
heraetian  conjugate  of  Equation  (C.21)  : 


(C.21) 

take  the 


(at  “  AJ*  *  0axatax**To 


(C.22) 


and  Multiply  on  the  right  by  If  we  choose  ta  such  that 


Va 


Vo 


(C.23) 


we  obtain 


<dt  -  c\)*  *  -  DVtax*0 


(C.24) 


Note  that  the  contrasting  foras  of  Equations  (C.21)  and  (C.24)  depend  on  the  use 
of  hypercoaplex  nuabers  and  the  relation  (C.23).  A  suitable  choice  for  rA  is 
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r. 


This  choice  yields  the  relations 


and 


Va 


Vk 


■H 
A  0 


-  V 


Va 


(c 

(•< 
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APPENDIX  D 


THE  SCALAR  STRESS  TENSOR 


Since  the  Lagrangian  density  contains  second-order  derivatives  of  the 
fields,  it  is  necessary  to  extend  the  standard  derivation  of  the  stress-energy 
tensors.  Denoting  general  derivatives  of  the  fields  f  and  9  by  and  the 
variation  in  the  Lagrangian  density,  4/f,  is  given  by 


Cl 


(D.  1 


In  the  following  discussion,  specific  Mention  of  the  adjoint  teres  will  be 
dropped,  i.e.,  the  existence  of  parallel  adjoint  teres  will  be  implied.  From 
Appendix  A,  the  minimum  in  the  variation  in  satisfies  an  equation  similar  to 
that  of  (A. 12): 


q  ■  — m  -  .  -A.  JL* 


Op 

where  xy  is  one  of  the  components  of  the  four  vector 
x  *  (r,t) 

Substitution  of  Equation  (D.2)  in  Equation  (D.1)  gives 


If  the  transformations  of  the  fields  are  such  that 
♦(x)  -  ♦,(x)  *  t(x)  ♦  <5*(x)  , 

and 

♦(x)  -  *^{x)  »  ^(x)  ♦  «*p(x)  , 


(D.2 

(D.3 

(D.4 

(0.5] 

(D.6] 


then  the  total  variation  in  the  fields,  including  transformation  to  new 
coordinates  will  be  gi/en  by 


«_<Mx)  •  *'  (x)  -  *(x). 


(D. 


Vp(x>  "  *p(x' }  "  ♦ptxj* 


(0.8) 


Consequently, 


«T*(x)  -  *fr(x')  ♦  I  |*-«xy  • 


v  v 


(D.9) 


and 


3*. 


Wx)  ■  4VX’>  fI  ife\  • 


V  V 


(0.10) 


Substituting  Equations  (0.9)  and  (0.10)  into  Equation  (0.4)  yields 


(0.11) 


The  coefficient  of  8x,  in  Equation  (0.11)  is  the  tensor  T  given  in  Equation 

V  VU 

(3.10)  of  Part  I.  For  pure  translations  of  t  and  equals  zero  and 
Equations  (3.1)  to  (3.4)  give  Equations  (3.16)  and  (3.17). 


1  Vi  ■  Vi-  v-1  •  v  rw*FM<  r-w"  rv  *  a ■  n»r  *  v  wire  UXV  »  vnt\T»  un  V.H  U^(UV  >\n  VAlWin  *fl  vnnv 


APPENDIX  E 

In  the  special  case  where  the  approximation  of  constant  vorticity, 

Vx/7*  0 

may  be  made,  the  interaction  Hamiltonian,  can  be  written  as 


Jf  *Jl  +Jt  *Jf 
112  3, 

where 


Jf,  =  (VxA  -  7xA)*(7*yo*  ♦  bQ^)  * 


1  8c2 


Jfz  *  -k£(VxA7q  +  7qVxA)*(VW»  -  W)  , 
8c 


andj/3  is  the  same  as  Jf2  in  Section  2. 


E.1  Revised  Canonical  Transforms 


Me  first  note  that 


**►-*•*-  iIX-p5^[(P  ♦  q)(hjnp  -  cJCq> 


♦  (p  -  q)(b+C+  ♦  npCq)] 


^o*  *  *V*  ■  ToSS-tlP“m>  *  *  <P"cq  *  <“cqKPC« 

pa  q 


(E.l) 


(E.2) 


(E.3) 


(E.4) 


(E.5) 


Sffl* 


h 


m 

KK 


(P«*«  +  tn  Hn!;t  ♦  n,c„n 


(E.6) 


Therefore, 


•  ««  *  k>‘W  *  *W»q*,|q 

-  [('  *  *>%tk  *  q»„  KL_C_ 


♦ 


t%q  -  10  *  k>"<q*k 

[(0  *  k"W  -  0“{q 


]V-T 

'  q+k 


(E. 


Proceeding  in 


H 


2 


the  si 


»  vein,  m  find 
Jf2 


al).(2q  ♦  k) 
k 


(n 


q+k^q 


♦ 


n_C _ ♦  clc _ ) 

*  q+k  q  q+k 


(E. 


It  is  convenient  to  add  Equations  (E.7)  and  (E.8)  to  obtain  a  quantity  denoted 
by  h4: 


H4  »  H1  +  H2 


icy. 


2 


—  - kx(a. 

Wlp  qk  u  ✓S 


k 


q+k  q  q+k  q  q  q+k  q  q+k 


where  the  vector  functions  of  q  and  k,  A,  B,  C,  and  0,  are  given  by 


A(q,  k)  -  (2q  ♦  k)ttak  ♦  (q  ♦  k)W||q+k  +  qu^  , 


B(q,  k)  -  (2q  +  k)^  +  (q  +  k)^  -  qu^  . 


C(q,  k)  -  (2q  ♦  k)(*ak  +  qu^  -  (q  +  k)Ufq+k  . 


0(q,  k)  *  (2q  +  k)o»ak  -  qufq  -  (q  ♦  k)ufq+k 


The  required  "canonical"  transforms  needed  to  diagonalize  H.  are 

4 


cm  *  HQ  *  °^q  *  "  q),7a  +  CQ  *  «  *  «  q)C+]//2U  , 

q  q  ak  nq  q  q  ak  Cq  q  q 

Ci  "  “V  “V  *  Vmi  *  ‘V  “.kq  '  “cqq>V/V*Iq  • 

dq  ■  11 V  “qkq  -  V>q,  *  <V  “,kq  *  “{,qKi>^q  ' 


di  •  t(V  V  -  WK  *  <V  *V  *  %qKq  ’^q  ’ 

The  resulting  contribution  to  Ha  will  lead  to  a  Hamiltonian  of  the  form 

P 


Ivf 

H»  ■  ^“**k*k  -  -J=  to,(\  -  *il  pj.kcq  -  dq*kdq>'“. 


■  P“.k*k*k  -  ,fV*k  *  ‘I'l  ,neqk  -  "dqk1/“,!  ' 


E.2  Order  Parameter 


An  alternative  means  of  gaining  insight  into  the  behavior  of  the  flow  system 
is  to  examine  the  behavior  of  variables  analogous  to  the  model  of  the  lasing 
transition  studied  by  Haken51  56  and  others. E7”51*  in  this  formulation,  a 
variable  related  to  a^  and  a£  at  low  values  of  k  will  behave  much,  like  an  order 
parameter  in  the  Landau-Ginzberg  model  of  critical  behavior.515  518  We  begin 
with  a  truncated  Hamiltonian,  omitting  and  H3  for  the  time  being.  We  find 
the  following  equations  of  motion  for  the  normal  modes  of  H^: 


a_  =  -iu  .a_  +  fk(kxe)X[AijJ  .n  +  Bn  +  Cn  C _  +  DflC _ ]  (E.20) 

k  aK  k  K  q  q+*  q  q  K  q  q  q+k  q  q+k 


aj  *  i«akat  +  fk(kxe)X[ArfJ+khq  +  Brjq+kC^  +  CnqC _  +  Dele _ ]  »  (E.21) 


q+k  q  q+k 


",  ’  ‘  l  V°*'\  *  ‘l'  '%-k  *  BC— 1  • 


(E.22) 


q+k 


nj  *  i»_  _n*  -  I  fkkx{a.  -  a+)[An*  k  +  cc _ ]  , 

q  nq  q  k  k  k  k  q+k  q+k 


(E . 23 ) 


*  l  v°«ak  -  *  °y  • 


(E.24) 


and 


£  *  1uc,cJ  *  l  V'x'N  -  *i)[cVk  *  DCk-,J 


(E.25) 


Equations  (E.22)  to  (E.25)  are  now  used  to  find  the  equations  of  motion  for 


N  and  N.  ,  as  well  as  for  two  new  variables,  A  and  a  ,  defined  by 

Hd  Cq  q  q 

<7  ■  n_C_ 

q  q  q 


(E.26) 
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Mot®  that 


CVV3  *  Vq.q'  ' 


(E.28 


tv*;.  ' 


(E.29 


(E.30 


The  equation  of  notion  for  N  Is  obtained  by  multiplying  Equation  (E.22)  by  p* 
and  Equation  (E.23)  by  and  adding.  We  adopt  the  convention  that  the 
appearance  of  terns  with  k  In  their  subscripts  Implies  a  sun  of  terms  with  both 
±|k| .  We  find 


n„  *  I  fkkx(ak  -  a!).[Cn  C _  -  • 


(E.31 


Similarly,  letting 


\  •  vx<%  -  <>■ 


(E.32 


■  l  'k-[CC-Vk  -  ■ 


(E.33 


|  VW{-V*  *  *  "J*V  *  %C-1 


(E.34 


,ccwq  -  l  V'**^  *  C(C5C^  *  "5Vk>  *  S?k-q! 


(E.3S 


We  now  take  the  1i.nit  k-o  and  find  with 


"nq  +  °Cq 


(E.36) 


the  following  approximate  equations: 


(E.37) 


*  -icqcrq  ♦  Jig  Fk.[(A+0)oq  ♦  Bnq] 


(E.38) 


and 


•5 


*  *V[(A+0)<TJ  +  Cnq3 


(E.39) 


Equations  (E.37)  to  (E.39)  can  be  used  as  the  starting  point  for  two 
different  treatments  of  critical  behavior.  The  first  is  a  computer  study  of  a 


transition  to  chaos  similar  to  that  of  the  Lorentz  model  of  Benard  convection. 


E19 


The  second  parallels  a  model  of  the  lasing  transition  studied  extensively  by 
E1-E6 


Haken . 
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APPENDIX  F 


K 


TRANSFORMATION  OF  THE  INTERACTION  HAMILTONIAN 

In  this  appendix,  we  show  how  Haailtonians  such  as  H.  given  in  Equation 

P 

(3.1)  nay  be  transformed  to  the  form  given  in  Equation  (3.6).  The  treatment 
follows  closely  that  given  by  Kittel  for  the  electron-phonon  interaction.F1 
Consider  the  truncated  Hamiltonian  H^,  given  by  a  diagonal  basis  Hamiltonian 
Hq  and  a  perturbation  part  H': 


where  for  example 


Vk\ 


(F.1) 


(F.2) 


and 

H'  -  to«.{  -  •  (F.3) 

The  key  step  in  transforming  Equation  (F.3)  is  to  postulate  the  existence  of  a 
function  S  such  that 


S  -  [S,Hd]  «  -  jp-  H'  .  (F.4) 

q 

Equation  (F.4)  can  be  shown  to  be  satisfied  using  the  interaction  picture. F1 

Although  we  are  discussing  a  fully  classical  system,  the  concept  of  Schroedinger, 

F2 

Heisenberg,  and  interaction  picture  can  still  be  applied  ,  as  implied  by  the 
discussion  in  Appendix  B. 


Integrating  Equation  (F4)  gives 


s«»  •  ~  -T~  /°dt  M'(t)  . 

q  -D 


(F.5) 


[134] 


Substituting  Equation  (F. 5)  into  Equation  (F.4)  and  then  substituting  the 
resulting  font  of  Equation  (F.4)  into  Equation  (F.1)  gives 


ht  •  HD  -  S1  £H’<o)'S(o)]  ' 

q 


(F.6) 


Fron  Equation  3.1 


fl2  Q 

H_  ■  K_  +  —  /  dt[H*(t),H’(o)] 
2u*  -0 

q 


f'  -  Vk  • 


(F.7) 


(F.8) 


Therefore,  the  perturbative  tern  in  Equation  (3.1)  in  the  fields  Cq  can  be  put 


in  the  forn  (y  ■  -1) 
o 


f2 

J°  dt[kx[ak(t)  -  a!(t)]cj+k(t)cq(t),k-x(ak,  -  )cj,+k,cq(  ].(F.9) 


2«;  -d 

q 


This  nay  be  written  as 


h=  -  z dt 

The  Fourier  transfornation  yields,  upon  contracting  the  Poisson  brackets, 

,  1  fkf k  \2  / ca*kcaCQ' -kco*  ca»kCQCa ' -kca '  \ 

C  “  2  \Uq  )  \  Wcq+k  -  %  +  wak  ucq+k  "  %  “  “ak  /  * 


(F . 10) 


(P.11 ) 


(135] 


S 


Performing  the  addition  and  making  a  similar  transformation  for  the  d  fields 
gives  Equation  (3.6). 

One  difficulty  occurs  in  that  there  are  several  nonlinear  contributions  to  H 
of  the  form  described  above.  Transformations  using  functions  S  of  the  type 
described  in  Equations  (F.4)  and  (F.S)  must  be  applied  to  the  various  nonlinear 
contributions  simultaneously.  This  can  be  achieved  through  the  use  of  projec¬ 
tion  operators,  again  using  the  concepts  described  in  Appendix  B. 

REFERENCE 
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1967). 


APPENDIX  Q 


In  this  appendix.  Me  will  display  the  steps  required  to  generate  the 
specific  tera  shown  in  Equation  (4.21)  of  Part  III.  We  begin  with  the  set  of 
Equations  (4.7)  to  (4.12).  For  the  particular  tera  of  interest,  we  replace  the 
0%  tera  in  Equation  (4.7)  with  expressions  (4.9)  and  (4.10).  We  find  that 

X  e^^A^q.kJn+fq+kWq)  * 

X  £n|i|k1A|i(q,k)Qn(q^OQn+(q)Fn(q+k)F,,(q) 
q 

-  X  I  £ili|k1AiB(q,k)6n(q+k)Gn+(q)Fn(q+k)f(k")erstks 
q  k" 

[at(k")  -  aJ(k")][Ar(q,k")n(q+k")  ♦  Br(q,k")ft(q+k") ] 

♦II  .k /  (q,k)Gn+(q)Gn(q+k)Ff,{ q)f(k’)6  k 
q  Ha  1  a  npq  p 

[aq(k* )  -  aq(k' ) 1 [  An(q+k,k' )q(q-*k-k' )  ♦  Bn(q+k,k* )Cf(q+k+k* )  ] 

♦  0(  f2(k)  ).  (G. 1 ) 

In  this  expression  we  will  exaaine  the  second  tera  on  the  right  side,  looking  at 
the  te^  proportional  to  Ar.  Using  Equation  (4.9)  once  again,  but  this  time 
truncating  at  the  linear  tera,  we  find 

I  X  f(k")Gn(q^)G^+(q)G^(q^■k")F,,(q♦k)F^(q+k,,) 

q  k" 

«il.V.«')erst^  Vk">V«'k"> 

♦  the  same  tera  with  at(k")  -  a*(k").  (G.2) 


[137] 


He  trill  now  take  advantage  of  this  aquation,  and  tap toy  tha  relation 

<Fn(q+k)Fn(q+k")>  -  «(q+k-q-k") 

-  «(k-k")  (6.3) 

leaving  us  Mith,  after  the  sua  over  k"  is  carried  out, 

<  ^(k)k1«r#tks^  Gn(q*>G^<q)®n(q*>Vq'k)Vq'k>  ]at(k)  >  (G.4) 

replacing  at  with  a  lowest  order  contribution,  we  find  that  Equation  (4.7) 
becoaes 

<  «j(k)  >  <  [o*j(h)  ♦ 

®m j (k,€Mlaklerstks^  Qn(q^)Qn+(‘l)Qr,(^) 

f(k)A|>(q,k)f(k)Ar(q,k)Q*t(k)  ]p*(itl  >  (6.5] 

Thus,  we  write  to  lowest  order, 

<  a.(k)  >  •  <  6*  (k)  F?(k)  ♦  ... 

J  J  (6.6] 

**her*  ^**j(k)  is  now  given  by  Equation  (23).  Also,  it  is  seen  in  Equation  (4.5) 
that  the  sua  over  internal  wave- vectors  contributes  a  tera  of  the  fora 

Jmr{k)  "  I  A  <G'k>A  (q.lOG'T(q*IOGr?+(q)G,,(q+k)  . 
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